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ABSTRACT 

We present foreground-redueed CMB maps derived from the full Planck data set in both temperature and polarization. Compared to the corre¬ 
sponding Planck 2013 temperature sky maps, the total data volume is larger by a faetor of 3.2 for frequeneies between 30 and 70 GHz, and by 1.9 
for frequencies between 100 and 857 GHz. In addition, systematic errors in the forms of temperature-to-polarization leakage, analogue-to-digital 
conversion uneertainties, and very long time constant errors have been dramatically reduced, to the extent that the cosmological polarization signal 
may now be robustly recovered on angular scales I > 40. On the very largest scales, instrumental systematic residuals are still non-negligible 
compared to the expeeted cosmological signal, and modes with ^ < 20 are aeeordingly suppressed in the current polarization maps by high-pass 
filtering. As in 2013, four different CMB eomponent separation algorithms are applied to these observations, providing a measure of stability with 
respect to algorithmic and modelling choices. The resulting polarization maps have rms instrumental noise ranging between 0.21 and 0.27 pK aver¬ 
aged over 55' pixels, and between 4.5 and 6.1 pK averaged over 3.'4 pixels. The cosmological parameters derived from the analysis of temperature 
power spectra are in agreement at the Icr level with the Planck 2015 likelihood. Unresolved mismatches between the noise properties of the data 
and simulations prevent a satisfactory description of the higher-order statistical properties of the polarization maps. Thus, the primary applications 
of these polarization maps are those that do not require massive simulations for accurate estimation of uneertainties, for instance estimation of 
cross-spectra and cross-correlations, or stacking analyses. However, the amplitude of primordial non-Gaussianity is consistent with zero within 2cr 
for all local, equilateral, and orthogonal eonfigurations of the bispectrum, including for polarization E-modes. Moreover, exeellent agreement is 
found regarding the lensing E-mode power spectrum, both internally among the various component separation codes and with the best-fit Planck 
2015 ACDM model. 
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1. Introduction 

This paper, one of a set associated with the 2015 release of 
data from the Planck^ satellite, presents maps of the Cosmic 
Microwave Background (CMB) anisotropies derived from the 
full Planck data set, for a total of 50 months of observa¬ 
tions from the Low Frequency Instrument (LFI) and 29 for 
the High Frequency Instrument (HFI) (Planck Collaboration I 
2015). This analysis updates the temperature-only analysis of 
the first 15.5 months of Planck observations discussed in Planck 
Collaboration XII (2014), and presents the first CMB polariza¬ 
tion maps derived from Planck observations. 

Much of the Planck analysis effort since the 2013 data 
release has revolved around understanding and reducing in¬ 
strumental systematic uncertainties. As summarized in Planck 
Collaboration I (2015), this work has been highly successful, re¬ 
ducing the net power from systematic errors in the HFI CMB 
channels by almost two orders of magnitude on large angular 
scales (Planck Collaboration XI 2015). The main contributions 
to these improvements have come from improved temperature- 
to-polarization leakage modelling, reduced Analogue-to-Digital 
Conversion (ADC) errors, and improved modelling of Very 
Long Time Constants (VLTCs; Planck Collaboration VI 2015; 
Planck Collaboration VIII 2015). With these improvements, the 
Planck observations are now sufficiently free of instrumental ar¬ 
tifacts to allow a robust determination of the CMB polarization 
anisotropies on intermediate and small angular scales, covering 
multipoles i > 20. However, as described both in this paper and 
in Planck Collaboration VIII (2015); Planck Collaboration XI 
(2015), residual systematics are still not negligible compared to 
the CMB signal on the very largest scales {I < 20), and these 
modes are therefore removed by a high-pass filter from the cur¬ 
rent maps. The rate of progress is still excellent, though, and 
updated all-scale maps with low large-scale systematics are ex¬ 
pected to be released in the near future. 

For temperature, the most significant improvement in the 
Planck 2015 analysis pipeline is absolute calibration based on 
the orbital CMB dipole rather than the Solar dipole (Planck 
Collaboration I 2015). This change, combined with a better 
understanding of both the Planck beams and transfer func¬ 
tions (Planck Collaboration VI 2015; Planck Collaboration VIII 
2015), has reduced the uncertainties in absolute calibration to a 
few tenths of a percent, and the agreement between LFI, HFI, 
and the Wilkinson Microwave Anisotropy Probe (WMAP) has 
improved to the level of the uncertainties (Planck Collaboration 
12015). 

The component separation efforts of the Planck 2015 re¬ 
lease are summarized in three papers. The current paper is 
dedicated to CMB extraction, and presents the main Planck 
2015 CMB maps in both temperature and polarization. Planck 
Collaboration X (2015) addresses astrophysical component sep¬ 
aration as implemented by Commander (Eriksen et al. 2004, 
2008), a Bayesian component separation algorithm, and presents 
a global model of the microwave temperature sky ranging from 
408 MHz to 857 GHz, including detailed maps of synchrotron, 
free-free, spinning dust, thermal dust, and CO emission. In ad- 
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dition, a few minor components (line emission around 90 GHz 
and the thermal Sunyaev-Zeldovich effect near the Coma and 
Virgo clusters) are included in the model, as are instrumental pa¬ 
rameters in the form of calibration, bandpasses, monopoles, and 
dipoles. The corresponding polarization model includes only 
synchrotron and thermal dust emission. Planck Collaboration 
XXV (2015) presents a detailed analysis of the foregrounds be¬ 
low -100 GHz in both temperature and polarization. For detailed 
descriptions of the various foreground components relevant for 
microwave component separation, see either Sect. 2 of Planck 
Collaboration XII (2014) or Sect. 4 of Planck Collaboration X 
(2015). 

The foreground amplitude relative to CMB polarization is 
such that effective foreground suppression is required for al¬ 
most any cosmological analysis, but the optimal approach de¬ 
pends sensitively on the topic in question. For instance, because 
the fluctuation power of diffuse polarized foregrounds decays 
as a power-law in multipole moment i (Page et al. 2007; Gold 
et al. 2011; Planck Collaboration X 2015; Planck Collaboration 
Int. XXX 2014), it is of greater importance for high-/’ CMB 
power spectrum and likelihood estimation to minimize noise 
sensitivity and to marginalize over unresolved point sources, 
than it is to model diffuse foregrounds with high accuracy. In 
this case, it is more convenient to parametrize the residual fore¬ 
grounds in terms of power spectrum models, and to marginal¬ 
ize over these in terms of a few global parameters, than to 
marginalize over a large number of per-pixel foreground parame¬ 
ters. The Planck 2015 CMB likelihood therefore employs cross¬ 
spectra coupled to simple harmonic space foreground modelling 
(Planck Collaboration XI 2015), rather than the detailed fore¬ 
ground modelling described in this paper. 

Similarly, because of the low — but non-negligible — level 
of residual instrumental systematics on large angular scales in 
the Planck 2015 data, the low-/ likelihood also implements a 
special purpose cleaning algorithm, in terms of a simple tem¬ 
plate fit including only the cleanest 30, 70, and 353 GHz chan¬ 
nels (Planck Collaboration VI 2015; Planck Collaboration VIII 
2015; Planck Collaboration XI 2015). This approach is similar 
to that adopted by the 9-year WMAP likelihood (Bennett et al. 
2013), and allows for easy propagation of uncertainties from cor¬ 
related noise in terms of full pixel-pixel noise covariance matri¬ 
ces. 

Other applications, however, such as gravitational lens- 
ing and Integrated Sachs-Wolfe (ISW) reconstructions (Planck 
Collaboration XV 2015; Planck Collaboration XXI 2015), 
constraints on isotropy and statistics (Planck Collaboration 
XVI 2015), searches for primordial non-Gaussianity (Planck 
Collaboration XVII 2015), and constraints on global geometry 
and topological defects (Planck Collaboration XVIII 2015), re¬ 
quire actual CMB maps, and these are all based on the products 
described in this paper. 

As in 2013, we apply four complementary CMB com¬ 
ponent separation algorithms to the Planck 2015 sky maps. 
In alphabetical order, these are 1) Commander, a paramet¬ 
ric pixel-based Bayesian CMB Gibbs sampler (Eriksen et al. 
2004, 2008); 2) NILC, a needlet-based internal linear combi¬ 
nation method (Basak & Delabrouille 2012, 2013); 3) SEVEN, 
which implements linear template fitting based on internal tem¬ 
plates in pixel space (Leach et al. 2008; Femandez-Cobos et al. 
2012); and 4) SMICA, a semi-blind spectral-matching algorithm 
fully defined in harmonic space (Cardoso et al. 2008; Planck 
Collaboration XII 2014). These codes were all applied to CMB 
temperature reconstruction in Planck Collaboration XII (2014), 
and have now been extended to polarization, as described in 
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Appendices A-D. In addition, each algorithm is applied to sev¬ 
eral subsets of the full data set, including half-ring, half-mission, 
and yearly data splits (Planck Collaboration I 2015). Comparing 
the resulting maps, both between algorithms and data splits, pro¬ 
vides a good understanding of both instrumental and algorithmic 
uncertainties. 

The paper is organized as follows. In Sect. 2 we describe 
the Planck 2015 data selection and pre-processing. In Sect. 3 
we briefly review the component separation methods, deferring 
mathematical details to Appendices A-D. In Sect. 4 we present 
the derived CMB maps. In Sect. 5 we quantify the residual fore¬ 
ground emission present in the maps by cross-correlation with 
foreground templates. In Sect. 6 we present angular power spec¬ 
tra and corresponding cosmological parameters. In Sects. 7 and 
8 we consider higher-order statistics and gravitational lensing. In 
Sect. 9 we summarize the main features and limitations of these 
maps, and provide recommendations on their applications. We 
conclude in Sect. 10. 

2. Data selection and pre-processing 

In this paper we use the full-mission Planck data (Planck 
Collaboration VI2015; Planck Collaboration VIII2015) and ac¬ 
companying simulations (Planck Collaboration XII 2015). The 
CMB temperature maps are derived using all nine frequency 
channels, from 30 to 857 GHz. One of our component separa¬ 
tion methods. Commander, additionally uses the 9-year WMAP 
temperature sky maps (Bennett et al. 2013) and a 408 MHz sur¬ 
vey map (Haslam et al. 1982). The CMB polarization maps are 
derived from the seven frequency channels sensitive to polariza¬ 
tion, from 30 to 353 GHz. In most cases, the component separa¬ 
tion methods use frequency channel maps as input, with the ex¬ 
ception of the temperature analysis performed by Commander, 
which uses maps from subsets of detectors in each frequency 
channel, as specified in Table 1 of Planck Collaboration X 
(2015). 

The primary foreground-reduced CMB maps are derived 
from full-mission maps, maximizing signal-to-noise ratio and 
minimizing destriping errors. In addition, a number of data splits 
are analysed to enable internal consistency checks and to make 
estimates of the properties of the CMB maps. These data splits 
can also be used in analyses where more than one map is re¬ 
quired as input. For the purposes of component separation, each 
subset of the data must contain the same combination of fre¬ 
quency channels as the full-mission data set. The following data 
splits have been analysed: 

“ maps from the first and second half of each pointing period 
(“half ring”; HRl and HR2); 

“ maps from odd and even years, consisting of year l-r3 maps 
for the LFI channels plus year 1 maps for the HFI channels, 
and year 2-r4 maps for the LFI channels plus year 2 maps for 
the HFI channels (YRl and YR2); 

- maps from the first and second half of the mission, consisting 
of year l-r2 maps for the LFI channels plus half-mission 1 
maps for the HFI channels, and year 3-1-4 maps for the LFI 
channels plus half-mission 2 for the HFI channels (HMl and 
HM2). 

The HFI maps in the half-ring data split have spurious corre¬ 
lations between them (Planck Collaboration VIII 2015). If these 
maps are used to estimate the noise level in the power spectrum, 
the correlations cause the estimate to be biased low. This was al¬ 
ready seen in the analysis of the nominal-mission CMB maps in 


the 2013 release (Planck Collaboration XII 2014), and the same 
effect is seen in the full-mission maps. The odd- and even-year 
data split does not include all of the data; it omits HFI survey 5. 
As a consequence, the HFI maps have more missing pixels than 
in the half-mission data split. For these reasons, the half-mission 
maps have been used as the primary data split for assessing the 
properties of the CMB maps and for further analysis. 

2.1. Data 

The frequency maps are described in detail in Planck 
Collaboration II (2015) for LFI and Planck Collaboration VII 
(2015) for HFI. The maps have several features that are rele¬ 
vant to component separation. We summarize them here, and 
the reader is recommended to consult the references for further 
details: 

- Monopole and dipole contributions from the CMB, CIB, and 
other astrophysical components are estimated and removed 
during mapmaking (Planck Collaboration VI 2015; Planck 
Collaboration VIII 2015). This has an effect on component 
separation, and each method treats monopoles and dipoles in 
a different way, as described in Appendices A-D. 

- The HFI maps are corrected for zodiacal light emission 
(ZLE) by subtracting a model of the emission at the ring level 
during mapmaking (Planck Collaboration VIII 2015). This 
differs from the treatment in the 2013 component separa¬ 
tion (Planck Collaboration XII 2014), where the ZLE model 
was not subtracted. 

- Leakage from intensity to polarization due to bandpass mis¬ 
matches between detectors is estimated based on the ground 
measurements of the bandpasses (Planck Collaboration IV 
2015; Planck Collaboration VI 2015; Planck Collaboration 
VIII 2015). These estimates are subtracted from the polar¬ 
ized maps. 

- The HFI maps at 100, 143, and 217 GHz are renormalized 
in order to correct for far-sidelobe effects in the calibra¬ 
tion (Planck Collaboration VII 2015). 

- Missing pixels are filled with the average values of pixels 
in the surrounding area. This area is defined as being within 
a radius of 1° for LFI maps and HFI frequency maps, and 
within a radius of L5 for HFI detector subset maps. 

Point source catalogues and masks have been provided as 
input to the component separation. Construction of the PCCS2 
catalogue is described in Planck Collaboration XXVI (2015). 
Masks have been constructed based on these catalogues. For 
each frequency channel, the intensity source mask removes all 
sources detected down to a S/N threshold of 4 for the LFI chan¬ 
nels and 5 for the HFI channels. For each polarized frequency 
channel, the polarization source mask removes all sources de¬ 
tected at 99.99 % significance or greater. The details of the mask¬ 
ing procedures used by each component separation method are 
described in Appendices A-D. 

2.2. Sinnulations 

To validate our results, we analyse realistic simulations of the 
Planck data set called full focal plane 8 (FFP8). They are based 
on detailed models of the instrument and sky, and are described 
in full in Planck Collaboration XII (2015). We summarize their 
contents here. 
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2.2.1. CMB 

The CMB was simulated using an input ACDM model based on 
the 2013 cosmological parameter results (Planck Collaboration 
XVI 2014). The fiducial simulation contains no primordial ten¬ 
sor modes or primordial non-Gaussianity. However, four variants 
of the same CMB realization have been produced that include 
non-zero values of the tensor-to-scalar ratio and non-Gaussianity 
of a local type. 

2.2.2. Foregrounds 

The Planck sky model (PSM) has been used to simulate the 
foreground components. The intensity part of the simulation in¬ 
cludes all astrophysical components that were identified in the 
2013 release. The diffuse components that are relevant at low fre¬ 
quencies consist of synchrotron, free-free, and anomalous dust. 
At high frequencies, CO, thermal dust, and CIB are included. 
The foreground modelling has been improved in this version of 
the simulations. In particular, Planck 353 GHz data were used to 
improve the frequency scaling of the dust emission. In polariza¬ 
tion, the diffuse foregrounds are synchrotron and thermal dust. 
The extragalactic emission from radio and infrared sources has 
been simulated in intensity and polarization, and the SZ effect 
from clusters of galaxies has been included in intensity. 

2.2.3. Simulated observations 

Time-ordered data (TOD) for each detector were simulated using 
the satellite pointing, and the individual detector beams, band- 
passes, noise properties, and data flags. The same mapmaking 
used for the data is used to generate maps from the simulated 
TOD. All of the maps from subsets of the data have also been 
generated from the simulations. 

Two versions of the maps are available, with and without 
bandpass mismatch leakage. The latter is simulated using the 
average bandpass for all detectors in a frequency channel, elimi¬ 
nating the leakage effect. The version of the maps without band¬ 
pass leakage is considered in this paper. 

In addition to the fiducial maps, a set of 10 000 Monte Carlo 
realizations of CMB and noise has been generated. These real¬ 
izations are intended to be used to assess the uncertainties on the 
results. 

2.2.4. Mismatch between simulations and data 

In analysing the simulations, a number of deficiencies have be¬ 
come evident. First, the amplitude of the CMB component does 
not match that of the data. This is the (expected) consequence of 
the fact that the CMB model for the simulations was specified 
before the recalibration of the Planck data between the 2013 re¬ 
lease and the present one was completed. This mismatch can 
be mitigated by increasing the amplitude of the CMB simula¬ 
tions by 1.3 % when comparing them to the data. Second, the 
noise properties of the simulated maps do not precisely match 
those of the data. This does not significantly affect the analysis 
of the CMB temperature maps, since they are signal-dominated. 
However, it does affect the analysis of CMB polarization maps, 
because they are more noisy. The noise mismatch appears to be 
scale-dependent, since the adjustment of the amplitude of the 
noise simulations to match the data depends on the resolution of 
the maps. This is explored in Sect. 7. 

Any analysis that relies on simulations to estimate the uncer¬ 
tainties of a result from the CMB polarization will be limited by 


these mismatches. Despite this, many analyses are possible, in¬ 
cluding those using cross-spectrum, cross-correlation, or stack¬ 
ing techniques. 

3. Component separation methods 

The four methods used by Planck to separate the CMB from 
diffuse foreground emission were described in detail in Planck 
Collaboration XII (2014). They are representative of the main 
approaches to component separation developed in recent years. 
The methods can be divided into two types. The first type 
assumes only knowledge of the blackbody spectrum of the 
CMB, and the foregrounds are removed by combining the multi¬ 
frequency data to minimize the variance of the CMB component. 
The second type constructs an explicit parameterized model of 
the CMB and foregrounds with an associated likelihood, and 
the CMB component is obtained by maximizing or sampling 
from the posterior distribution of the parameters. Either type of 
method may be implemented in the map domain or in the har¬ 
monic domain. We recall briefly their main features and com¬ 
ment on their application to polarization data. Descriptions of 
the changes in each algorithm since 2013 are given in the appen¬ 
dices. 

- Commander (Eriksen et al. 2006, 2008) is a Bayesian para¬ 
metric method that works in the map domain. Both the CMB 
and foregrounds are modelled using a physical parameter¬ 
ization in terms of amplitudes and frequency spectra, so 
the method is well suited to perform astrophysical com¬ 
ponent separation in addition to CMB extraction (Planck 
Collaboration X 2015). The joint solution for all compo¬ 
nents is obtained by sampling from the posterior distribu¬ 
tion of the parameters given the likelihood and a set of pri¬ 
ors. To produce a high-resolution CMB map, the separation 
is performed at multiple resolutions with different combi¬ 
nations of input channels. The final CMB map is obtained 
by combining these solutions in the spherical harmonic do¬ 
main. This obviates the need for the Ruler step that was 
used in 2013 to extend the Commander solution to high res¬ 
olution. A low-resolution version of the separation is used 
to construct the temperature power spectrum likelihood for 
large angular scales, as described in Planck Collaboration XI 
(2015). Note that Commander employs detector and detector 
set maps rather than full frequency maps, and excludes some 
specific detector maps judged to have significant systematic 
errors, in addition to incorporating the 9-year WMAP tem¬ 
perature sky maps and a 408 MHz survey map. Thus, the se¬ 
lection of data is not identical between Commander and the 
other three methods. 

- NILC (Delabrouille et al. 2009a) is an implementation of in¬ 
ternal linear combination (ILC) that works in the needlet 
(wavelet) domain. The input maps are decomposed into 
needlets at a number of different angular scales. The ILC so¬ 
lution for the CMB is produced by minimizing the variance 
at each scale. This has the advantage that the weights used 
to combine the data can vary with position on the sky and 
also with angular scale. The solutions are then combined to 
produce the final CMB map. 

- SEVEM (Fernandez-Cobos et al. 2012) is an implementation 
of the template cleaning approach to component separation 
that works in the map domain. Foreground templates are typ¬ 
ically constructed by differencing pairs of maps from the 
low- and high-frequency channels. The differencing is done 
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Table 1. Masks and statistics of component-separated CMB maps from data and FFP8 simulations. 






Method 


Parameter 


Coimnander 

NILC 

SEVEM 

SMICA 

Sky fraction, /s^y [%]“ 

Data confidence mask 

T. . . . 

81.9 

96.4 

84.5 

85.0 


Q,u . 

83.1 

96.5 

79.4 

85.0 

preferred mask 

FFP8 confidence mask 

T . . . . 
Q.U.. 

T . . . . 

75.3 

96.9 

77.6 

77.4 

82.8 

86.5 


Q,U.. 

87.5 

96.1 

79.3 

83.4 

preferred mask T ... . 

Q.U .. 

Standard deviation at FWHM 10', Nside 

Data confidence mask T ... . 

= 1024 [jjK]'’ 

101.8 (4.4) 

101.6 (5.4) 

73.5 

75.7 

101.4 (3.2) 

101.1 (4.2) 


e.... 

6.3 (5.8) 

5.3 (4.8) 

6.3 (6.3) 

5.2 (4.7) 


u .... 

6.3 (5.8) 

5.3 (4.7) 

6.3 (5.9) 

5.2 (4.5) 

preferred mask 

T . . . . 

101.3 (4.4) 

100.9 (5.3) 

101.3 (3.2) 

101.0 (4.1) 


e .... 

6.3 (5.8) 

5.2 (4.7) 

6.3 (6.3) 

5.2 (4.7) 


U ... . 

6.3 (5.8) 

5.2 (4.5) 

6.3 (5.9) 

5.2 (4.5) 

FFP8 confidence mask 

T . . . . 

104.5 (3.5) 

106.5 (4.6) 

104.0 (4.5) 

104.3 (3.5) 


e .... 

5.6 (5.0) 

5.1 (4.4) 

6.1 (5.6) 

5.0 (4.2) 


U ... . 

5.7 (5.1) 

5.2 (4.4) 

6.1 (5.6) 

5.0 (4.3) 

preferred mask 

T ... . 

104.3 (3.5) 

107.5 (4.2) 

104.3 (4.4) 

104.5 (3.5) 


Q . 

5.6 (5.0) 

5.0 (4.3) 

6.1 (5.6) 

4.9 (4.2) 


u ... . 

5.7 (5.1) 

5.0 (4.3) 

6.1 (5.6) 

5.0 (4.2) 

Standard deviation at FWHM 160', Nside = 64 

Data confidence mask T.... 48.0 (1.30) 

48.7 (1.00) 

47.5 (0.37) 

47.5 (0.79) 


e.... 

0.25 (0.21) 

0.31 (0.27) 

0.29 (0.26) 

0.29 (0.26) 


u . . . . 

0.25 (0.21) 

0.31 (0.26) 

0.29 (0.25) 

0.28 (0.25) 

preferred mask 

T . . . . 

47.4 (1.29) 

47.1 (1.01) 

47.3 (0.37) 

47.1 (0.78) 


e .... 

0.25 (0.21) 

0.30 (0.27) 

0.29 (0.27) 

0.28 (0.26) 


U ... . 

0.25 (0.21) 

0.30 (0.26) 

0.29 (0.25) 

0.28 (0.25) 

FFP8 confidence mask 

T . . . . 

55.3 (0.34) 

59.8 (0.83) 

55.4 (0.48) 

55.6 (0.50) 


e .... 

0.23 (0.18) 

0.27 (0.22) 

0.27 (0.22) 

0.24 (0.19) 


U ... . 

0.23 (0.18) 

0.27 (0.23) 

0.26 (0.23) 

0.34 (0.20) 

preferred mask 

T ... . 

55.4 (0.35) 

61.7 (0.79) 

55.6 (0.48) 

56.0 (0.49) 


Q . 

0.23 (0.18) 

0.26 (0.22) 

0.27 (0.22) 

0.24 (0.19) 


u ... . 

0.23 (0.18) 

0.26 (0.22) 

0.26 (0.23) 

0.24 (0.20) 


^ Sky fractions are given at A^side = 2048 for T and A^side = 1024 for Q and U. 

^ Values in brackets are standard deviations of half-mission half-difference maps, giving an indication of the level of residual noise and systematic 
effects. Standard deviations of Q and U have been computed from high-pass filtered maps. For details of the downgrading and high-pass filtering 
procedures, see text. 


in order to null the CMB contribution to the templates. These 
templates are then used to clean each CMB-dominated fre¬ 
quency channel by finding a set of coefficients to minimize 
the variance of the map outside of a mask. Thus SEVEN pro¬ 
duces multiple foreground-cleaned frequency channel maps. 
The final CMB map is produced by combining a number of 
the cleaned maps in harmonic space. 

- SMICA (Cardoso et al. 2008) is a non-parametric method that 
works in the spherical harmonic domain. Foregrounds are 
modelled as a small number of templates with arbitrary fre¬ 
quency spectra, arbitrary power spectra and arbitrary cor¬ 
relation between the components. The solution is obtained 
by minimizing the mismatch of the model to the auto- and 
cross-power spectra of the frequency channel maps. From 
the solution, a set of weights is derived to combine the fre¬ 


quency maps in the spherical harmonic domain to produce 
the final CMB map. Maps of the total foreground emission 
in each frequency channel can also be produced. In the anal¬ 
ysis performed for the 2013 release (Planck Collaboration 
XII2014), SMICA was the method that performed best on the 
simulated temperature data. 

3.1. Extension to polarization 

The methods described above were applied to Planck tempera¬ 
ture data for the 2013 release (Planck Collaboration XII 2014), 
and they have been extended to operate on polarization data in 
the present work. A key distinction between the methods is the 
choice of operating domain. Two of the methods, Coimnander 
and SEVEM, operate in the map domain, so it is most natural for 
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UT78 


UPB77 



Fig. 1. Preferred masks for analysing component-separated CMB maps in temperature {left) and polarization {right). 


them to do the polarized component separation on the Q and U 
maps. The other two methods, NILC and SMICA, operate in the 
harmonic or needlet domain. An intrinsic part of the transform of 
polarized maps to these domains is the decomposition of Q and 
U into E and B modes, which is accomplished by using spheri¬ 
cal harmonic transforms on the full sky. Thus these two methods 
perform their separation directly on E and B. 

3.2. Outputs 

In addition to producing CMB maps, each method provides 
“confidence” masks to define the region of the sky in which the 
CMB solution is trusted in temperature and polarization. The 
procedure each method uses to define the mask is described 
in Appendices A-D. The confidence masks are used to define 
masks for further analysis of the data. Two of the pipelines, 
Coinmander and SMICA also produce foreground products, de¬ 
scribed in Planck Collaboration X (2015). 

The first 1000 Monte Carlo realizations of CMB and noise 
have been propagated through the four pipelines. This has been 
done twice, once using the parameters derived from the data, and 
once using the parameters derived from the fiducial FFP8 maps, 
to provide a set of simulations to accompany each data set. 

The methods produce maps at different resolutions, as de¬ 
scribed in Appendices A-D. The products have been brought 
to a standard resolution to compare them and for distribution. 
Standard resolution temperature maps have a Gaussian beam of 
5' FWHM and HEALPix^ resolution Aside = 2048. Standard res¬ 
olution polarization maps have a Gaussian beam of 10' FWHM 
and HEALPix resolution Aside = 1024. If maps are produced at 
higher resolution then they are downgraded to these standard res¬ 
olutions. The downgrading procedure for maps is to decompose 
them into spherical harmonics (T, or £” and B, as appropriate) 
on the full sky at the input HEALPix resolution. The spherical 
harmonic coefficients, < 2 ^^, are convolved to the new resolution 
using 


a 


out 

tm 




,out 






( 1 ) 


lution is first downgraded like a temperature map. The smooth 
downgraded mask is then thresholded by setting pixels where 
the value is less than 0.9 to zero and all others to unity to make 
a binary mask. This has the effect of enlarging the mask to ac¬ 
count for the smoothing of the signal. In addition to the standard 
resolution products, the maps, masks, and Monte Carlo realiza¬ 
tions have been downgraded to lower resolutions for analyses 
that need them, using the above procedure. 

The polarization maps and Monte Carlo realizations have 
been decomposed into E and B mode maps, and downgraded to 
lower resolutions too, for analyses that work on E and B directly. 
The decomposition is done on the full-sky maps using spherical 
harmonic transforms. The CMB maps from the real-space meth¬ 
ods, Commander and SEVEM, are inpainted before doing the de¬ 
composition. Both, the standard Commander and SEVEM CMB 
maps are inpainted inside their corresponding confidence masks 
using a constrained realization (see Appendix A for details). For 
both methods, only the CMB maps are inpainted, not the Monte 
Carlo realizations as this would be too computationally expen¬ 
sive. The other two methods, NILC and SMICA, work on E and 
B modes, so it is possible to make E and B maps directly from 
their outputs in addition the the standard Q and U maps. 


4. CMB maps 

In this section we present and discuss the component separated 
CMB maps in temperature and polarization. For temperature, we 
compare the maps to those extracted from the nominal mission 
data in 2013. We compare the maps from different methods to as¬ 
sess their consistency, and we use the FFP8 simulations to assess 
the accuracy of the methods and the robustness of the solutions. 
Throughout the discussion, we make use of appropriate masks in 
order to highlight differences at intermediate and high Galactic 
latitudes, ignoring the plane of the Galaxy where differences are 
much higher due to the complexity of the foreground signal and 
its dominance over the CMB. 

4.1. Temperature maps 


where is the beam transfer function, is the HEALPix pixel 
window function, and the “in” and “out” superscripts denote the 
input and output resolutions. They are then synthesized into a 
map directly at the output HEALPix resolution. Masks are down¬ 
graded in a similar way. The binary mask at the starting reso- 

^ http://healpix.sourceforge.net 


Temperature confidence masks produced by the methods have 
been used to make combined masks for further analysis of 
the data. The first mask is constructed as the union of the 
Commander, SEVEM, and SMICA confidence masks. The NILC 
mask is not included in the union because it removes a signif¬ 
icantly smaller fraction of the sky. This union mask has /sky = 
11.6%. We refer to it as UT78 and adopt it as the preferred mask 
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Fig. 2. Component-separated CMB temperature maps at full resolution, FWHM 5', = 2048. 

Commander 2013-2015 NILC 2013-2015 



SEVEM 2013-2015 SMICA 2013-2015 
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Fig. 3. Differences between the component-separated CMB temperature maps from the 2013 and the 2015 releases. The maps have 
been smoothed to FWHM 80' and downgrading to NsHe = 128. 
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Fig. 4. Pairwise difference maps between CMB temperature maps. As in the previous Fig. 3, the maps have been smoothed to 
FWHM 80' and downgraded to Aside = 128. 
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Fig. 5. Difference between output and input CMB temperature maps from FFP8 simulations. Smoothing and downgrading as in 
Fig. 4. 
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for analysing the temperature maps. It is shown on the left in 
Fig. 1. An extended version of the mask has been constructed 
by adding to the UT78 mask those pixels where the the standard 
deviation between the temperature maps is greater than lOyuK. 
This mask has /sky = 76.1 %, and we refer to it as UTA76. A 
union mask has been created for the FFP8 simulations in the 
same way as for the data. It has /sky = 73.5 %, and we refer to it 
as FFP8-UT74. 

The CMB temperature maps produced by the four methods 
are shown in Fig. 2. No obvious differences are seen in these 
maps, and they appear visually consistent outside the mask. An 
important assessment of the robustness and consistency of the 
CMB T component separation solutions is provided in Fig. 3, 
which shows the differences between the Planck 2013 and 2015 
maps for each method. Several interesting features may be seen 
in these differences, most of which correspond directly to a 
better understanding of the systematic uncertainties in the new 
maps. Starting with Commander, the most striking features are 
large-scale swaths tracing the Planck scanning strategy with a 
peak-to-peak amplitude of ~ lOyuK. This pattern is very simi¬ 
lar to that originally pointed out by Larson et al. (2014), who 
found this by subtracting the 9-year WMAP ILC map (Bennett 
et al. 2013) from a template-cleaned version of the Planck 2013 
100 GHz map. Similar patterns are also seen in the Commander 
residual maps shown in Fig. 2 of Planck Collaboration X (2015), 
corresponding to detector maps that are rejected from the new 
2015 analysis. These structures are primarily due to two ef¬ 
fects, namely destriping errors from bandpass mismatch between 
detectors and far sidelobe contamination (Planck Collaboration 
VIII 2015; Planck Collaboration X 2015). By rejecting partic¬ 
ularly susceptible channels in the updated Commander analysis, 
these errors are greatly reduced in the new map. 

Turning to the other three difference maps in Fig. 3, we 
see that the residuals are internally very similar, but quite dif¬ 
ferent from the Commander residuals. In these cases, the two 
most striking features are, first, a ~ 5 yuK quadrupole roughly 
aligned with the CMB solar dipole (Galactic coordinates (/, b) = 
(264°,48°); Planck Collaboration I 2015), and, second, clear 
traces of the zodiacal light emission. The former is explained 
by the fact that the HFI data processing in 2013 did not subtract 
the relativistic Doppler quadrupole of ~ 6 yuK, which is now sub¬ 
tracted in the 2015 maps. Similarly, the latter is explained by the 
fact that ZLE was not subtracted in the mapmaking in 2013, but 
is subtracted in the updated processing (Planck Collaboration VI 
2014; Planck Collaboration VIII 2015). Commander, on the other 
hand, is less sensitive to residual ZLE, for the following two rea¬ 
sons. First, in 2013 it used channels only up to 353 GHz, which 
are less affected by the ZLE than the higher frequencies. Second, 
by virtue of fitting independent thermal dust spectral parameters 
(index and dust temperature) per pixel, it can efficiently absorb 
the ZLE in the thermal dust component. However, some of the 
ZLE may still be observed in the Commander differences: rem¬ 
nants of the “red arcs” typical in the second and fourth quadrant 
of the sky are just visible in the difference of the Commander 
solutions, while being very evident in all other cases. 

The residuals seen in Fig. 3 are small compared to the typi¬ 
cal CMB anisotropies, with features mostly constrained to < 5- 
10 yuK, with a distinct large-scale pattern. In particular, these 
small differences are completely negligible for power spectrum 
and cosmological parameter estimation. The only cosmological 
application for which some care is needed is the study of large 
scale isotropy (Planck Collaboration XXIII 2014), for instance 
with respect to CMB quadrupole-octopole alignment. An up¬ 
dated isotropy analysis of the new sky maps will be presented in 


Planck Collaboration XVI (2015); no significant differences are 
found compared to the 2013 results. 

Table 1 summarizes main properties of the CMB maps de¬ 
rived both from the data and from the fiducial set of FFP8 sim¬ 
ulations. We evaluate standard deviations in two cases, corre¬ 
sponding to high (FWHM 10', Aside = 1024) and intermediate 
resolution (FWHM 160', Aside = 64). The values in parentheses 
are standard deviations of half-mission half-difference (HMHD) 
maps, and they give an estimate of the level of uncertainties due 
to instrumental noise and systematic effects. The same quantities 
are given for the FFP8 simulations. At this level, results show 
good consistency for both data and simulations. The SEVEN maps 
have the lowest standard deviation, as measured by the HMHD 
maps, at small and intermediate angular scales. 

Pairwise differences between all four maps are shown in 
Fig. 4, after smoothing to 80' FWHM and downgrading to 
Aside = 128. As expected, differences are largest close the edge 
of the mask, where the absolute foreground level is the highest. 
Comparing these with the corresponding maps from 2013 shown 
in Fig. 6 of Planck Collaboration XII (2014), and noting that the 
new color bar spans a range that is four times narrower than the 
previous one, we see that the internal agreement between the 
four methods is substantially better in the new maps, typically 
by about a factor of two. 

Figure 5 shows the differences between the FFP8 outputs and 
the input CMB map. The residuals are smallest for NILC and 
largest for SEVEN. However, we note that the foreground model 
adopted for the simulations was chosen to be more complex than 
the real sky, in order to explicitly probe modelling errors. In 
particular, the simulated thermal dust frequency spectrum ex¬ 
hibits a strong positive (and spatially dependent) curvature at 
low frequency that is neither captured in the parametric mod¬ 
els adopted by Commander, nor easily modelled by the spatial 
templates adopted by SEVEN. This additional complexity makes 
it hard to draw a strong conclusion about the performance on 
the real data, for which the thermal dust spectrum may be very 
well approximated by a simple one-component greybody com¬ 
ponent with a nearly spatially constant spectral index (Planck 
Collaboration X 2015). 

4 . 2 . Polarization maps 

We now turn our attention to the foreground-reduced 
CMB polarization maps. As discussed extensively in Planck 
Collaboration I (2015), Planck Collaboration VI (2015), and 
Planck Collaboration VIII (2015), the residual systematics in the 
Planck polarization maps have been dramatically reduced com¬ 
pared to 2013, by as much as two orders of magnitude on large 
angular scales. Nevertheless, on angular scales greater than 10°, 
correponding ioi < 20, systematics are still non-negligible com¬ 
pared to the expected cosmological signal. Different combina¬ 
tions of input frequency channels for the component separation 
have been explored in order to mitigate the polarization residu¬ 
als. However, it was not possible, for this data release, to fully 
characterize the large-scale residuals from the data or from sim¬ 
ulations. Therefore the CMB polarization maps provided in the 
current release have been high-pass filtered to remove the large 
angular scales. This has been implemented by applying a cosine 
filter to the E and B spherical harmonic coefficients of the maps. 
This filter is defined as 





qi-cos(;r|E|) 

1 <i <^2 5 

(2) 

U 
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Fig. 6. Component-separated CMB Q maps at resolution FWHM 10', A^side = 1024. 
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Fig. 7. Component-separated CMB U maps at resolution FWHM 10', A^side = 1024. 
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Fig. 8. 20° X 20° patch of the high-pass filtered Commander CMB polarization map, centered on the North Ecliptic Pole, (/, b) = 
(96°, 30°). Each map is pixelized with a HEALPix resolution of A^side = 1024, and has an angular resolution of 10' EWHM. The top 
row shows Q and U maps derived from the full-mission data set, the middle row shows the corresponding E and B maps, and the 
bottom row shows the E and B maps of the half-ring half-difference (HRHD) map. Note the characteristic -h and x patterns in the 
Q and U maps, and the clear asymmetry between E and B in the full data set. Also note that the HRHD E map is consistent with 
both the full and HRHD B maps. 
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Fig. 9. Pairwise differences between CMB Q maps, after smoothing to FWHM 80' and downgrading to A^side = 128. 


and we have used = 20 and I 2 = 40. The same filtering has 
been applied to the FFP8 fiducial maps and to the Monte Carlo 
simulations. 

As for temperature, individual polarization confidence 
masks derived for each method have been used to make com¬ 
bined masks for further analysis. Our first polarization mask is 
simply the union of the Commander, SEVEM, and SMICA confi¬ 
dence masks, which has /sky = 77.6% and we refer to it as 
UP78. However, the 1-point statistics analysis summarized in 
Sect. 7.1, revealed significant point source contamination in the 
Commander, NILC, and SMICA maps using this mask. SEVEM was 
not affected by this problem because it applies an inpainting 
technique to remove the brightest point sources (see Appendix C 
for further details). For this reason, two extended versions of the 
mask were created, the first by excluding in addition the pixels 
where the standard deviation between the CMB maps, averaged 


in Q and U, exceeds 4yuK. This mask has /sky = 76.7 %, and 
we refer to it as UPA77. The second was made by additionally 
excluding from the union mask the polarized point sources de¬ 
tected at each frequency channel. It has /sky = 77.4 %, and we 
refer to it as UPB77. This mask is shown on the right of side 
of Fig. 1, and we adopt this as the preferred polarization mask, 
since it is physically better motivated than UPA77. Also, although 
it keeps a larger fraction of the sky than UPA77, it is sufficient to 
alleviate the point source contamination. 

Masks have been made in a similar way for the FFP8 sim¬ 
ulations. The union mask has /sky = 76.3 %, and we refer to 
as FFP8-UP76. An extended mask that also excludes polar¬ 
ized point sources has /sky = 75.7 %, and we refer to it as 
FFP8-UPA76. FFP8-UPA76 is the preferred mask for FFP8 po¬ 
larization analysis. 
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Fig. 10. Pairwise differences between CMB U maps, after smoothing and downgrading as in Fig. 9. 


Figures 6 and 7 show the high-pass filtered Q and U Stokes 
parameters of the CMB maps after applying the UPB77 mask. 
The maps are shown at full resolution, and are thus dominated 
by instrumental noise except in the regions at the ecliptic poles 
where integration time is greatest. Visually, the methods operat¬ 
ing in the harmonic (NILC and SMICA) and spatial (Commander 
and SEVEM) domains are more similar to each other than the 
other methods. 

In order to have better visual insight at the map level, in 
Figure 8 we show a 20° x 20° patch of the high-pass-filtered 
Commander polarization maps centred on the North ecliptic pole. 
The top row shows the full-mission Q and U maps. Note the 
characteristic “-f” pattern in Q and “x” pattern in U ; this is the 
expected signal for a pure E mode signal. To make this point 
more explicit, the middle row shows the same map decomposed 
into E and B components. There is a clear asymmetry between 


them, with E having visibly more coherent power than B, again 
as expected for an ^’-dominated signal. Finally, the third row 
shows the half-ring half-difference (HRHD) E and B maps, il¬ 
lustrating the noise level in the full-mission maps. Comparing 
the middle and bottom rows, there is clearly an £'-mode excess 
in the full-mission map, whereas the corresponding full-mission 
5-mode map is consistent with the HRHD 5-mode map. In ad¬ 
dition, the HRHD E’-mode map is also consistent with both the 
full-mission and HRHD 5-mode maps, suggesting that all are 
consistent with instrumental noise. 

Pairwise differences between the four polarization maps are 
shown in Figs. 9 and 10. The NILC and SMICA solutions appear 
closest to each other. The regions of the sky that are most af¬ 
fected by differences appear to be those with a higher noise level, 
as may be seen by comparing to Figs. 6 and 7. For completeness. 
Figs. 11 and 12 show the differences between the FFP8 outputs 
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Fig. 11. Difference between output and input CMB Q maps from FFP8 simulations. Smoothing and downgrading as in Figs. 9 and 
10 . 
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Fig. 12. Difference between output and input CMB U maps from FFP8 simulations. Smoothing and downgrading as in Fig. 11. 
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Table 2. Correlation coefficients of CMB maps with foreground templates for temperature and polarization. 


Correlation Coefficient 


Foreground Template 

Commander 

NILC 

SEVEN 

SMICA 

TEMPERATURE 

Ha . 

0.010 + 0.071 

0.011+0.071 

0.019 + 0.071 

0.003 + 0.057 

CO. 

-0.004 + 0.027 

-0.003 + 0.027 

-0.003 + 0.027 

-0.007 + 0.022 

857 GHz. 

-0.043 + 0.084 

-0.032 + 0.084 

-0.037 + 0.084 

-0.029 + 0.083 

Haslam. 

-0.062 + 0.115 

-0.051 + 0.116 

-0.065 + 0.115 

-0.023 + 0.069 

POLARIZATION 

WMAP K-Ka. 

-0.057 + 0.026 

-0.116 + 0.024 

-0.026 + 0.025 

-0.027 + 0.026 

WMAP K-Ka (HPF) . . . 

0.0042 + 0.0036 

0.0054 + 0.0037 

0.0147 + 0.0037 

0.0092 + 0.0036 


and the input CMB map. The differences show a pattern similar 
to that of the noise, though with higher amplitude with respect 
to pairwise differences of solutions from data, possibly reflecting 
again the enhanced complexity of the simulated sky with respect 
to real data. 

The combination of high-pass Altering and noise makes vi¬ 
sual comparison of these maps difficult. The rms summary pro¬ 
vided in Table 1 is more informative in this respect. Comparing 
the HMHD rms values listed in parentheses, we see that NILC 
and SMICA have the lowest effective polarization noise levels at 
high angular resolution, with rms values that are roughly 20 % 
lower than those observed for Commander and SEVEN. One plau¬ 
sible explanation for this difference is the angular resolution 
adopted for the fitting process by the four methods; whereas 
Commander and SEVEN perform the polarization analysis at 10' 
FWHM resolution, NILC and SNICA adopt a 5' FWHM resolu¬ 
tion. When comparing the rms values at an angular resolution 
of 10', as presented in Table 1, the maps from the latter two 
methods are smoothed by post-processing to a lower resolution, 
whereas the maps from the two former codes are not. 

This effect is not relevant at intermediate angular scales, for 
instance at 160' FWHM, as shown in the bottom section of 
Table 1. On these angular scales, we see that the situation among 
the codes is reversed, and Commander provides a 20 % lower ef¬ 
fective noise than the other three codes. 


5. Correlation with externai tempiates 

Correlation of the CMB maps with templates of foreground 
emission provides a first diagnostic of residual contamination in 
the maps. We compute the correlation coefficient r between two 
maps X and y as 


1 y X,- - <x) ji - <y) 

^pix “ 1 ■ tTx (Ty 


(3) 


where the index i runs over the A^pix pixels observed with the 
common mask, (x) = 2, x,VA?pix, ctx = [Z,(Xi - (x))^/(AVix - 
1)]*^^, and similarly for y. Both maps and templates are 
smoothed to FWHM 1° and downgraded to A^side = 256 before 
computing the correlation, and the analysis is performed sepa¬ 
rately on temperature and polarization maps. 

The foreground templates considered for temperature are: 
the 408 MHz radio survey of Haslam et al. (1982); the velocity- 
integrated CO map of Dame et al. (2001); the full-sky Ha tem¬ 
plate of Finkbeiner (2003); and the Planck 857 GHz channel 


map. The uncertainty in the value of r due to chance correla¬ 
tions between foregrounds and the cleaned maps is estimated by 
computing the correlations between the templates and the 1000 
simulations of CMB and noise provided by each method. For po¬ 
larization, the only template we consider is one constructed from 
the WMAP 9-year maps. The template is made by smoothing the 
K and Ka band maps to FWHM 1° and differencing them to re¬ 
move the CMB contribution. This produces a template contain¬ 
ing the low-frequency polarized foreground emission and noise. 
The correlation analysis is done twice, once with the original 
maps and templates, and once with a high-pass Altered version. 
The resulting coefficient factors and 1 a uncertainties are shown 
in Table 2. 

For temperature, the results for all methods are compati¬ 
ble with no correlation within Icr. From this, we conclude that 
there are no significant temperature foreground residuals with 
the same morphology as the templates in the map, to a precision 
set by cosmic variance. For polarization, the analysis of unfll- 
tered maps show that SEVEN and SNICA are compatible with no 
correlation at the 1 cr level. Commander has a moderate level of 
correlation at the 2 a level, and for NILC we And a correlation 
around 4.5 cr. For high-pass Altered maps (labelled by HPF), the 
level of correlation is reduced for Commander and NILC to about 

1.5 cr, while it increases for SEVEN and SNICA to about 4 and 

2.5 cr, respectively. These detections may be due to accidental 
correlations of the map with the template which are not taken 
into account when computing the uncertainty. Neither the noise 
in the template nor the systematics in the maps at large angular 
scales are modelled, the latter being important for the unflltered 
case. In the Altered case, the signal is reduced relative to the 
noise in the template, which could exacerbate spurious correla¬ 
tions. From this analysis, it is difficult to draw strong conclusions 
about the residual contamination in the polarized CMB maps. 

6. Power spectra and cosmological parameters 

In this section we evaluate the foreground-cleaned maps in terms 
of CMB power spectra and cosmological parameters. We em¬ 
phasize that the Planck 2015 results parameters are not based on 
the high-resolution foreground-cleaned CMB maps presented in 
this paper, but are instead derived from the likelihood described 
in Planck Collaboration XIII (2015). That likelihood combines 
the low-resolution Commander temperature map derived from 
Planck, WMAP, and 408 MHz with a template-cleaned LFI 
70 GHz polarization map in a pixel-based low-^ likelihood, and 
adopts a cross-spectrum-based estimator for the high-^ temper- 
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ature and polarization likelihood. The high-^ likelihood relies 
on a careful choice of masks along with templates and mod¬ 
elling, ah in the power spectrum domain, to reduce the contribu¬ 
tion from diffuse Galactic and discrete Galactic and extragalac- 
tic foreground emission. The templates and source models are 
marginalized over when estimating cosmological parameters. 

The parameter estimation methodology that we use here is 
primarily a tool to evaluate the quality of the high-resolution 
CMB maps and to assess overall consistency with the Planck 
2015 likelihood. We start with the foreground-cleaned CMB 
maps and masks described in Sect. 4. The maps have been 
cleaned from diffuse foregrounds and, to a varying extent, from 
extragalactic foregrounds. We construct simplihed templates for 
the residual extragalactic foregrounds that are marginalized over 
when estimating cosmological parameters. 

While the Planck 2015 likelihood takes into account calibra¬ 
tion and beam uncertainties, we have not done so in this anal¬ 
ysis. Two of the methods, Coinmander and SMICA, ht for the 
relative calibration between frequency channels, but the uncer¬ 
tainties from this process are not propagated into the maps. The 
other two methods, NILC and SEVEN, assume that the frequency 
channels are pefectly calibrated. None of the four methods prop¬ 
agate the uncertainties in the beam transfer functions into the 
CMB maps. 

We are interested in assessing the relative quality of the CMB 
maps, for which it is more important to assess the spread of pa¬ 
rameters between methods and as a function of angular scale, 
rather than to provide absolute numerical values. However, since 
the CMB maps we describe here are the basis for the analysis of 
the statistical isotropy of the CMB, primordial non-Gaussianity, 
and gravitational lensing by large scale structure, it is both im¬ 
portant and reassuring that the parameter values that we hnd are 
reasonably close to those of Planck Collaboration XIII (2015). 

Figure 13 shows the TT and EE power spectra of the 
foreground-cleaned maps, applying the UT78 mask in temper¬ 
ature and the UP78 mask in polarization. The power spectra 
of the half-mission half-sum (HMHS) data are shown as thick 
lines; those of the half-mission half-difference (HMHD) data are 
shown as thin lines. The HMHS spectra contain signal and noise, 
whereas the HMHD spectra contain only noise and potential sys¬ 
tematic effects. The variations in the temperature noise spectra 
at low to intermediate multipoles are caused by the component 
separation methods optimizing the trade-off between foreground 
signal and noise as a function of scale. At high mutipoles, the 
same spectra are smooth because the relative weighting of the 
frequency channels is set by the noise levels. The breaks in 
the Commander noise spectra are caused by the hybridization of 
maps at different resolutions to make the hnal map. 

Figure 14 shows the power spectra of the signal in the maps 
estimated using XFaster (Rocha et al. 2010, 2011) from the 
HMHS and HMHD maps. The top panels compare each of the 
four power spectra derived from the component separation maps 
with the best-ht ACDM power spectrum derived from the Planck 
likelihood including multipoles up to ^ = 2500. The bottom pan¬ 
els show the spectrum differences between the component sepa¬ 
ration methods and the best-ht spectrum. The bottom left panel 
shows signihcant differences between the temperature spectra 
computed from the four raw component separation maps. The 
Commander spectrum has signihcantly more temperature power 
than the other three methods at ^ > 1000, with a clear break 
in amplitude between i = 1000 and 1200. As described in 
Appendix A, Commander only uses frequencies 217 GHz and 
above for multipoles above i ~ 1000; the Commander hybridiza¬ 
tion approach achieves greater angular resolution at the price of 


a higher point source contribution above i = 1000, compared to 
the other three methods. On the other hand, htting foregrounds 
pixel-by-pixel ensures that point sources remain localized in the 
hnal map, and the shape of the foreground spectrum has the 

cx shape expected for a Poisson source contribution. 

SMICA has the lowest high-^ temperature power excess from 
extragalactic sources of ah four methods. In this case, though, 
the excess is almost hat between i ^ 1200 and 1700, rehect- 
ing the effective harmonic space frequency weighting as a func¬ 
tion of multipole, and the effective source contribution there¬ 
fore has a more complicated overall behaviour. To account at 
least partially for this effect, we use the FFP8 simulations to 
construct an effective extragalactic source template as a func¬ 
tion of multipole for each method individually, and marginalize 
over the corresponding amplitude during parameter estimation. 
These templates, however, are only as good as the inputs used 
for the simulations, and uncertainties in the frequency depen¬ 
dency of the underlying true source populations translate into 
an ^dependent error in the source templates for NILC SEVEN, 
and SMICA. Commander is more robust against this particular ef¬ 
fect because it smooths ah maps to a common angular resolution 
prior to component separation, but, as seen in Fig. 14, this ro¬ 
bustness comes at the price of a higher overall source amplitude. 

The bottom right panel of Fig. 14 shows the corresponding 
information for the EE spectrum. Once again, we see a clear 
positive excess in these residual spectra. In the absence of a de¬ 
tailed FFP8 model for this excess, we assume for now that it 
has a shape oc and marginalize over its amplitude during 
cosmological parameter estimation. 

Cosmological parameters from the component separated 
maps in both temperature and polarization are determined using 
XFaster power spectra, coupled to CosmoMC (Lewis & Bridle 
2002) using a correlated Gaussian likelihood. Specihcahy, we 
include multipoles between ^Ynin “ ^max, where ^max = 

1000,1500, or 2000 for temperature, and 1000 or 1500 for polar¬ 
ization. We adopt a standard six-parameter ACDM model, and, 
since low-T” data are not used in the likelihood, impose an infor¬ 
mative Gaussian prior of r = 0.07 + 0.006. As mentioned above, 
we construct foreground templates for each CMB map by propa¬ 
gating the simulated full-sky FFP8 foreground maps through the 
respective pipeline and estimating the resulting power spectra 
normalized to some pivotal multipole. 

The resulting cosmological parameters are summarized in 
Fig. 15 for both TT (hhed symbols) and EE (unhhed symbols). 
Starting with the temperature cases, we hrst observe good over¬ 
all internal agreement between the four component separation 
methods, with almost ah differences smaller than ~ 1 cr within 
each multipole band. Second, we also observe good agreement 
with the best-ht Planck 2015 ACDM model derived from the 
likelihood, as most of the differences are within 1 cr. The no¬ 
table exception is the power spectrum amplitude, Ase~^^, which 
is systematically low at ~ 2 cr for ^max = 1000 for ah methods. 

On a more detailed level, however, there is some evidence of 
internal tensions at the 1-2 cr level, most clearly seen in 
and For these two parameters, there are almost 1 cr shifts 

for NILC, SEVEN, and SMICA going from ^max = 1000 to 1500, 
and again from ^max = 1500 to 2000. Commander appears to be 
somewhat more robust with respect to multipole range than the 
other three methods. Although some variation is indeed expected 
by statistical variation alone, the combination of the shapes of 
the power spectrum differences seen in Fig. 14 and the system¬ 
atic parameter trends suggest systematic uncertainties at the 1- 
2 cr level due to extragalactic foreground modelling, as discussed 
above. For this reason, we do not recommend using the compo- 
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Fig. 13. Power spectra of the the foreground-cleaned CMB maps. Left. TT power spectra evaluated using the UT78 mask. Right EE 
power spectra evaluated using the UP78 mask. The thick lines show the spectra of the half-mission half-sum maps containing signal 
and noise. The thin lines show the spectra of the half-mission half-difference maps, which given an estimate of the noise. The black 
line shows the Planck 2015 best-fit CMB spectrum for comparison. 
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Fig. 14. CMB TT (left) and EE {right) power spectra for each of the four foreground-cleaned maps. Top panels show raw bandpow- 
ers with no subtraction of extragalactic foregrounds; the grey lines show the best-fit ACDM model from the Planck 2015 likelihood. 
The bottom panels show residual bandpowers after subtracting the best-fit ACDM model showing the residual extragalactic fore¬ 
ground contribution. The foregrounds are modelled and marginalized over when estimating parameters, see Figure 15. 
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Fig. 15. Comparison of cosmological parameters estimated from the TT and EE spectra computed from the foreground-cleaned 
CMB maps. Within each group, the three left-most points show results for TT with = 1000, 1500, and 2000 the two right-most 
points show results for EE with = 1000 and 1500. For comparison, we also show the corresponding parameters obtained 
with the Planck 2015 likelihood including multipoles up to ^max = 2500 as the horizontal line surrounded by a grey band giving 
the uncertainties. The foreground model used for the cleaned CMB maps is the method-tailored full-sky model from the FFP8 
simulations. 


nent separated maps for cosmological parameter estimation at 
this time; for this purpose the Planck likelihood method is pre¬ 
ferred, which shows much better stability with respect to multi¬ 
pole range (Planck Collaboration XI 2015). 

For polarization, the results are more ambiguous, with fluc¬ 
tuations relative to the temperature prediction beyond 2 cr. Note, 
however, that the cosmic variance contributions to the temper¬ 
ature and polarization parameters are essentially independent, 
and the two estimates are therefore only expected to agree sta¬ 
tistically, not point-by-point. Still, in particular SEVEM appears 
to show evidence of larger deviations than expected in polariza¬ 
tion for several parameters, and the 6 parameter for Commander 
in polarization shows some hints of tension with respect to the 
corresponding temperature estimate. 

Overall, however, we And good consistency between the four 
different component separation methods in both temperature and 
polarization. For temperature, the main outstanding issues are 
uncertainties in the residual extragalactic foreground model at 
the ~ lOyuK^ level for NILC, SEVEM, and SMICA, and a substan¬ 
tially larger effective point source amplitude in the Commander 
map compared to the other three maps. The parameters derived 
from polarization observations are generally in good agreement 
with the corresponding temperature parameters, except for the 
outliers noted above. 


7. Higher-order statistics 

We now consider the higher-order statistics of the CMB maps in 
the form of 1-point statistics (variance, skewness, and kurtosis), 
A^-point correlation functions, and primordial non-Gaussianity 
We focus in particular on the polarization maps and their 
degree of consistency with the FFP8 simulations. The tempera¬ 
ture results are described in Planck Collaboration XVI (2015), 
Planck Collaboration XVII (2015), and Planck Collaboration 
XVIII (2015). 

7.1. I-point statistics 

The variance, skewness, and kurtosis of the maps, and the pre¬ 
processing steps needed to compute them, are described in detail 
in Planck Collaboration XVI (2015), Monteserm et al. (2008), 
and Cruz et al. (2011). The procedure is to normalize the data, 
dp, in pixel p by its expected rms, dp = dpIcTp. The rms is cal¬ 
culated from 1000 FFP8 realizations (Planck Collaboration XI 
2015), for both CMB anisotropies and instrumental noise. To 
the extent that represents an accurate description of the data 
variance, and both the sky signal and the instrumental noise are 
Gaussian-distributed quantities, d will be Gaussian distributed 
with zero mean and unit variance. In temperature, the variance 
of the instrumental noise is subtracted in order to determine the 
variance of the CMB signal in the data. Once the variance of the 
CMB signal is estimated, it is used to extract the skewness and 
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kurtosis from the normalized map. This procedure is well es¬ 
tablished, and provides a direct test for the presence of residual 
foregrounds and of CMB Gaussianity. 

The temperature analysis reveals that the Planck 2015 CMB 
maps and the FFP8 fiducial CMB maps are fully compatible 
with the Monte Carlo simulations, therefore they can be used 
for further statistical analyses. For more details about the tem¬ 
perature results and the FFP8 validation analysis please refer to 
Appendix E.2 and Planck Collaboration XVI (2015). 

For polarization, since the Q and U maps are not rotationally 
invariant, we consider the polarized intensity P = ^JQ^ + P 
is not Gaussian-distributed with zero mean, and its skewness and 
kurtosis are non-vanishing; however, by comparing the data with 
the Monte Carlo ensemble, we can test for residual foregrounds 
and quantify the performance of the component separation meth¬ 
ods. Noise in the P maps is complicated; rather than trying to 
remove the noise contribution (as is done in the case of tempera¬ 
ture), we compare the P maps with the Monte Carlo simulations 
of CMB and noise together. Table 3 gives the resulting lower- 
tail probabilities, that is, the percentage of simulations that show 
a lower variance, skewness, or kurtosis than the P maps. This 
is done at three different resolutions (A^side = 1024, 256, and 64) 
for both the data maps (using mask UPB77) and the fiducial FFP8 
maps (using mask FFP8-UPA76). 

At lower resolutions. Table 3 shows good agreement be¬ 
tween the fiducial CMB maps and the Monte Carlo simulations 
for all methods. At high resolution, skewness, and kurtosis are 
mostly in agreement with the simulations, with NILC and SMICA 
showing slightly high values for the kurtosis (lower tail probabil¬ 
ities of 98.5 % and 99.4 %, respectively). However, there is ex¬ 
cess variance in the maps at high resolution, 2-3.5 a away from 
the mean of the simulations, with Commander deviating the least. 

Using the individual components of the FFP8 fiducial maps 
(CMB, noise, and foregrounds), we are able to quantify the con¬ 
tributions to the statistics of the high-resolution maps from each 
component separately. Fig. 16 compares the variance (left col¬ 
umn), skewness (middle column), and kurtosis (right column) 
extracted from the FFP8 Monte Carlo simulation P maps (his¬ 
tograms) with those extracted from the FFP8 fiducial realiza¬ 
tion of: the sum of CMB and noise (blue); the sum of CMB, 
noise, and thermal dust (green); the sum of CMB, noise, and ra¬ 
dio point sources (orange); and the sum of CMB, noise, and all 
foregrounds (red). The last defines the lower tail probabilities in 
Table 3. We have also investigated other components, but found 
that their contributions are negligible. 

We see that the only case that is fully compatible with the 
Monte Carlo ensemble for all methods is that of CMB and noise 
(blue lines). In particular, adding the thermal dust component 
(green lines) increases the variance slightly outside the accept¬ 
able range for both NILC and SMICA and the same effect be¬ 
comes even stronger when adding the rest of the foreground 
components for these methods. Coimnander is only slightly af¬ 
fected by foregrounds by this measure, and remains within the 
acceptable range even for the full foreground model, while 
SEVEM is an intermediate case. For the kurtosis, we find rel¬ 
atively high sensitivity to radio source residuals (the orange 
and red lines coincide), but very low sensitivity to diffuse fore¬ 
grounds. We conclude that the anomalous statistics seen in anal¬ 
ysis of the high-resolution component-separated FFP8 CMB 
maps are due to the foreground components, which could plau¬ 
sibly be caused by the additional complexity of the FFP8 fore¬ 
ground model with respect to the real sky. The Monte Carlo sim- 


Table 3. Percentage of simulations showing a lower variance, 
skewness, or kurtosis than the P maps, the “lower-tail probabil¬ 
ity”. The top half of the table shows results for the data, and 
the bottom half for the FFP8 fiducial map, using the UPB77 and 
FFP8-UPA76 masks, respectively. In brackets, in the variance 
column, we report the excess variance of the data in percent with 
respect to the mean of the variance distribution of the Monte 
Carlo simulations. 



Lower-tail Probability 

Map 

Variance 

Skewness 

Kurtosis 

Data 

A^side = 1024 

Commander .. 

100.0 (3.5) 

89.5 

45.9 

NILC. 

100.0 (4.0) 

99.9 

95.9 

SEVEM . 

100.0 (4.5) 

93.6 

85.1 

SMICA . 

100.0 (3.2) 

73.0 

49.7 

A^side - 256 

Commander .. 

100.0 (7.5) 

62.3 

44.5 

NILC. 

100.0 (8.6) 

54.8 

39.4 

SEVEM . 

100.0 (9.2) 

57.6 

52.8 

SMICA . 

100.0 (6.9) 

33.0 

24.2 

A^side — 64 

Commander .. 

100.0(11.7) 

0.0 

0.1 

NILC. 

100.0 (18.0) 

22.3 

22.1 

SEVEM . 

100.0 (20.5) 

14.6 

6.6 

SMICA . 

100.0 (16.5) 

26.4 

30.3 

FFP8 

A^side — 1024 

Commander .. 

95.9 

77.0 

88.1 

NILC. 

100.0 

90.1 

98.5 

SEVEM . 

100.0 

64.6 

65.6 

SMICA . 

100.0 

91.2 

99.4 

A^side - 256 

Commander .. 

92.3 

72.5 

65.8 

NILC. 

82.5 

25.7 

23.1 

SEVEM . 

79.9 

17.1 

25.6 

SMICA . 

62.0 

62.3 

58.9 

A^side “ 64 

Commander .. 

97.5 

89.8 

73.0 

NILC. 

90.9 

53.0 

27.5 

SEVEM . 

68.7 

48.2 

33.1 

SMICA . 

78.1 

66.8 

35.9 


ulations are compatible with the CMB and noise components of 
the fiducial map, as they should be by design. 

Figure 17 shows the variance of the Monte Carlo simula¬ 
tions (from which the lower tail probabilities are summarized in 
the top half of Table 3) compared with that of the Planck 2015 
CMB P maps. The maps analysed are the sum of CMB and noise 
signal. Since the CMB is practically the same for all methods, 
the different average values of the variance of the Monte Carlo 
simulations are a reflection of somewhat higher noise on small 
angular scales in the SEVEM and Commander maps than in the 
SMICA and NILC maps. Looking at the data (red bars), we im¬ 
mediately see that they do not match the simulations. No simu¬ 
lation has a variance as high as the data for any method at any 
resolution. As already mentioned, this discrepancy is due to an 
underestimation of the noise in the FFP8 simulations (Planck 
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Variance Skewness Kurtosis 


Fig. 16. Polarised intensity variance (left column), skewness 
(middle column), and kurtosis (right column) evaluated from 
the FFP8 Monte Carlo simulations (histogram) and from com¬ 
ponents of the fiducial FFP8 map at A^side = 1024 outside the 
FFP8-UPA76 mask. The variance distributions have been nor¬ 
malized to the mean value of the Monte Carlo distributions for 
visualization purposes. Coloured vertical lines correspond to dif¬ 
ferent combinations of components: the sum of CMB and noise 
is shown in blue; the sum of CMB, noise, and thermal dust is 
shown in green; the sum of CMB, noise, and radio point sources 
is shown in orange; the sum of CMB, noise, and all foregrounds 
is shown in red. 



12 14 16 18 0.58 0.62 0.66 0.70 0.025 0.030 0.035 0.040 

Variance{fj,K^) 


Fig. 17. Polarized intensity variance evaluated from the FFP8 
Monte Carlo simulations (histogram) and from the Planck 2015 
maps (vertical red lines) outside the UPB77 mask. Columns from 
left to right show different resolutions (A^side = 1024, 256, and 
64), while rows show results for the four component separation 
methods. Unlike in Fig. 16, the variance distributions are not 
normalized. 


Collaboration XII 2015). In the variance column of Table 3, in 
parentheses, we report the excess variance of the data in percent 
relative to the mean variance of the Monte Carlo simulations. 
The excess is 3-4% at A^side = 1024, increasing to 10-20% at 


A^side = 64. Moreover, we see that there are some extreme values 
for both the skewness and kurtosis, in particular at low resolu¬ 
tion for Coininander and at high resolution for NILC. However, as 
long as the variance is inaccurate, it is difficult to decide whether 
residual foregrounds, a true non-Gaussian feature in the map, or 
noise underestimation is the cause. For now, we simply conclude 
that the simulations are inconsistent with the data, and this ef¬ 
fectively prevents studies of higher-order statistics of this kind 
for the polarization maps. For temperature, the same is not true 
because of the much higher signal-to-noise ratio, although care 
is warranted even then when probing into the noise-dominated 
regime above i ^ 1500-2000. 

We also performed the 1-point analysis on the Planck 2015 
polarization maps using the UP7 8 mask. In this case, the values 
of the skewness and kurtosis for SMICA, NILC, and Commander 
were significantly affected by the presence of point sources. 
Conversely, the results for SEVEN were consistent with the ex¬ 
pected distribution, because this method applies an inpainting 
technique to remove the signal from the brightest point sources 
(see Appendix C for more details). This motivated the construc¬ 
tion of the UPB77 mask, that excluded the brightest point sources 
detected in polarization, which strongly alleviated this problem 
for the other three methods as can be seen in the results of 
Table 3. The same behaviour was also found for the FFP8 fidu¬ 
cial maps. 

7.2. N-point correlation functions 

Real-space A/^-point correlation functions are a useful diagnostic 
of the statistics of CMB maps complementary to harmonic anal¬ 
yses. In this section we describe their application to the Planck 
2015 CMB polarization maps. Results for the FFP8 CMB maps 
are given in Section E.3. Details of their application to tempera¬ 
ture maps may be found in Planck Collaboration XVI (2015). 

For observed fields X measured in a fixed relative orientation 
on the sky, the V-point correlation function is defined as 


C^(^1,..., ^2v-3) = {X(hi ). • • X(h^)), (4) 

where the unit vectors hi,...,hj^ span an V-point polygon on 
the sky. Assuming statistical isotropy, A-point functions are 
functions only of the geometrical configuration of the A-point 
polygon. In the case of the CMB, the fields X correspond to AT 
and two Stokes parameters Q and U describing the linear po¬ 
larization of the radiation in direction h. In standard CMB con¬ 
ventions, Q and U are defined with respect to the local meridian 
of the spherical coordinate system of choice. However, Q and U 
form a spin-2 field and depend on a rotational coordinate sys¬ 
tem transformation. To obtain coordinate-system-independent 
A-point correlation functions the Stokes parameters are rotated 
with respect to a local coordinate system defined by the centre 
of mass of the polygon (see Gjerlpw et al. 2010 for details). The 
Stokes parameters in this new “radial” system are denoted by Qr 
and Ur^ 

Given rotationally invariant quantities X e {/ST, Qr, Ur], the 
correlation functions are estimated by simple product averages 
over all sets of A pixels fulfilling the geometric requirements set 
by 6 i,..., 62 N- 3 , which characterize the shape and size of the 
polygon. 


Cn(0i, ..., O2N-3) ■ 






(5) 
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The pixel weights Wj, • • •, are introduced to reduce noise or 
mask boundary effects. Masks set weights to 1 for included pix¬ 
els and 0 for excluded pixels. 

The shapes of the polygons selected for the analysis are 
pseudo-collapsed and equilateral configurations for the 3-point 
function, and a rhombic configuration for the 4-point function 
comprising two equilateral triangles sharing a common side. The 
4-point function is only computed in the analysis of tempera¬ 
ture maps. We use the same definition of pseudo-collapsed as 
in Eriksen et al. (2005), that is, an isosceles triangle where the 
length of the baseline falls within the second bin of the sepa¬ 
ration angles. The length of the longer edges of the triangle, 0, 
parametrizes its size. Analogously, in the case of the equilateral 
triangle and rhombus, the size of the polygon is parametrized by 
the length of the edge, 0. Note that these functions are chosen 
because of ease of implementation, not because they are better 
suited for testing Gaussianity than other configurations. In the 
following, all results refer to the connected 4-point function. 

We analyse the high-pass filtered Planck 2015 CMB maps 
at resolution FWHM 160', Aside = 64. We used the downgraded 
version of the UP78 mask in the analysis. 

The A-point functions are used to test the quality of the CMB 
estimates derived from the Planck data, and are shown in Fig. 18. 
We show the differences between the A-point functions for the 
high-pass filtered CMB maps and the corresponding mean val¬ 
ues estimated from the 1000 Monte Carlo simulations. The prob¬ 
abilities of obtaining values of the;^^ statistic for the Planck fidu¬ 
cial ACDM model at least as large as for the CMB map are given 
in Table 4. The results of the analysis of the temperature maps 
can be found in Planck Collaboration XVI (2015). 

The results for the data deviate significantly from the Monte 
Carlo simulations for almost all A-point functions involving at 
least two polarization fields. The smallest deviation is seen for 
the Commander map. The A-point functions for Monte Carlo 
simulations have smaller variance than for data. By compar¬ 
ing HMHD maps with Monte Carlo simulations of noise cor¬ 
responding to the CMB estimates, we find that the amplitude of 
the noise is underestimated by 18 %. After adjusting the Monte 
Carlo simulations to compensate, the A-point functions are more 
consistent with data. It is difficult to draw conclusions about the 
Gaussianity of the polarization maps from these results, since the 
A-point functions themselves are used to estimate the mismatch 
between the noise Monte Carlo simulations and the data. 

7.3. Primordial non-Gaussianity 

Primordial non-Gaussianity is often measured in terms of the 
amplitude, of the quadratic corrections to the gravita¬ 

tional potential, as well as by means of the three-point cor¬ 
relation function based on different triangle configurations. 
The results from these calculations for the foreground-cleaned 
CMB maps are presented in Planck Collaboration XVII (2015). 
Compared to the previous data release, we can now include both 
temperature and polarization bispectra in the analysis. We thus 
consider bispectrum estimates obtained from temperature 
and polarization data only, as well as the full constraint from all 
eight possible TTT, TTE, EET, and EEE combinations. 

Results obtained from application of the Komatsu-Spergel- 
Wandelt (KSW) estimator to the CMB maps after subtraction 
of the lensing-ISW correlation (see below) are listed in Table 5 
for various geometrical configurations that have been consid¬ 
ered (Planck Collaboration XVII 2015). It is interesting to eval¬ 
uate the impact of polarization data on the estimation of /nl 


Table 4. Probability-to-exceed (PTE) in percent for the A-point 
correlation function statistic applied to the Planck 2015 maps 
at Aside = 64 for each of the four methods, as shown in Fig. 18. 



PTE [%] 


Function 

Commander NILC SEVEM 

SMICA 


Two-point 


TQr . 

39.8 

7.1 

11.0 

6.4 

TUr . 

65.9 

39.4 

0.5 

31.5 

QrQr . 

<0.1 

0.4 

<0.1 

<0.1 

QrUr . 

<0.1 

0.4 

<0.1 

<0.1 

UrUr . 

<0.1 

0.4 

<0.1 

<0.1 

Pseudo-collapsed three-point 




TTQr . 

64.4 

61.2 

72.2 

72.1 

TTUr . 

44.2 

88.3 

69.6 

74.7 

TQrQr . 

21.9 

10.5 
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measurements. By considering only the temperature bispec¬ 
trum for the SMICA map, we obtain - 1-3 ± 5.7, 

while the polarization alone yields ± 31.0. 

Uncertainties are evaluated by means of Gaussian FFP8 simula¬ 
tions. We find consistency between pipelines at the 1 cr level. The 
results confirm and extend to polarization the absence of evi¬ 
dence of non-Gaussianity of primordial origin estimated through 
the bispectrum. The performance of the methods is also tested 
using Gaussian and non-Gaussian FFP8 simulations, showing 
that SMICA and SEVEM give the results closest to the inputs in 
both cases (see section 7.3 of Planck Collaboration XVII 2015). 

An interesting case to consider is that of the ISW-lensing 
bispectrum. Unlike with primordial shapes, there is a specific 
prediction for the expected amplitude of the ISW-lensing three- 
point signal in a given cosmological model. We can therefore 
use this shape to check whether the expected level of NG is re¬ 
covered, verifying that different component separation methods 
do not either spuriously add or remove any signal, at least in the 
squeezed limit where this shape is peaked. This is indeed the 
case. By normalizing the ISW-lensing shape in such a way as 
to have an amplitude of 1 in the best-fit model, we recover 
/nF(SMICA) “ 0.85 + 0.2 from the full analysis including all 
bispectra, /nF(SMICA) “ from temperature alone, and 

/nF(SMICA) “ from polarization alone. 

The constraint from polarization alone is, as expected, much 
looser than the one from temperature alone. However, polariza¬ 
tion does not have a negligible impact on the final combined 
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Fig. 18. The difference between the A^-point functions for the high-pass filtered A^side = 64 Planck 2015 CMB estimates and the 
corresponding means estimated from 1000 Monte Carlo simulations. The Stokes parameters Qr and Ur were locally rotated so that 
the correlation functions are independent of coordinate frame. The first row shows results for the 2-point function, from left to right, 
TQr, QrQr, ^nd QrUr- The second row shows results for the pseudo-collapsed 3-point function, from left to right, TTQr, TQrQr, 
QrQrUr, and UrUrUr, and the third row shows results for the equialteral 3-point function, from left to right, TTQr,TQrQr, QrQrUr 
and UrUrUr. The red solid, orange dot dot dot-dashed, green dashed and blue dot-dashed lines correspond to the Commander, 
NILC, SEVEN, and SMICA maps, respectively. The shaded dark and light grey regions indicate the 68% and 95% confidence regions, 
respectively, estimated using SMICA simulations. See Sect. 7.2 for the definition of the separation angle 6 . 


Table 5. Amplitude of primordial non-Gaussianity, /nl, es¬ 
timated by the KSW estimator. See Table 10 in Planck 
Collaboration XVII (2015) for full details. 


/nl 


Type Commander NILC SEVEN SMICA 


T local. 4+6 3+6 4+6 3+6 

T equilateral. -20+ 71 -28+ 69 -2+ 69 -11 + 70 

T orthogonal. -29+ 35 -45 + 33 -36+ 33 -34+ 33 

TE local. 4+5 1+5 -3+5 1+5 

TE equilateral. 14+ 46 -9+ 44 8+ 47 3+ 43 

TE orthogonal. -29+ 22 -25 + 21 -39+ 23 -25 + 21 

E local. 33 + 39 -1 + 33 60+ 42 26+ 32 

E equilateral. 327 + 165 75 + 140 292+ 167 144+ 141 

E orthogonal. -52+ 88 -78+ 76 -183 + 91 -128+ 72 


measurement, mostly due to contributions coming from TTE 
configurations. Adding polarization reduces the final uncertainty 
by ~ 30 %, as well as moving the recovered amplitude parame¬ 
ter closer to its expected value of 1. The implications of these 
results in terms of the physics of the early Universe, as well 
as the study of many additional shapes, are discussed in Planck 
Collaboration XVII (2015), which explains the algorithmic de¬ 


tails and all results from the procedure summarized here, and in 
Planck Collaboration XX (2015), which discusses the implica¬ 
tions for inflationary physics. 

8. Gravitational lensing by large-scale structure 

Gravitational lensing from intervening matter imprints a non- 
Gaussian signature in the CMB temperature and polarization 
maps, which in turn can be exploited to extract the gravitational 
potential integrated along the line of sight back to the surface 
of last scattering. Planck accurately measures the lensing poten¬ 
tial over most of the sky (Planck Collaboration XV 2015). As it 
remaps the CMB polarization, lensing partially transforms pri¬ 
mordial E modes into B modes, resulting in a secondary B mode 
spectrum peaking at around i = 1000. By forming a weighted 
product of the E modes from component separated CMB maps 
and the reconstructed lensing potential, it is possible to gener¬ 
ate a map of the expected lensing-induced B modes of CMB 
polarization. Cross-correlating this lensing 5-mode template to 
the total observed B modes provides an indirect measurement 
of the lensing 5-mode power spectrum (Planck Collaboration 
PIP116 2014; Planck Collaboration XV 2015). Lensing 5-mode 
template maps were synthesized for the four component separa¬ 
tion methods considered in this paper, using the Stokes parame¬ 
ter maps and the lensing potential reconstruction from the inten- 
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Fig. 19. Lensing-induced 5-mode power spectra in the 
component-separated polarization CMB maps. The solid line 
represents the best fit cosmology from the Planck data release 
in 2015. Error bars were evaluated using a semi-analytical ap¬ 
proximation validated over the FFP8 simulations as described in 
Planck Collaboration PIP116 (2014). 


sity map as discussed in Planck Collaboration PIP116 (2014). 
A common mask was generated by combining the union po¬ 
larization mask and the lensing potential 80 % mask of Planck 
Collaboration PIP116 (2014). After apodization using a cosine 
function over 3°, the resulting mask preserves an effective sky 
fraction of about 60 %. 

The lensing 5-mode power spectra for Coimnander, NILC, 
SEVEM, and SMICA, which have been measured by cross- 
correlating with the corresponding foreground cleaned polariza¬ 
tion maps, are shown in Fig. 19. The four CMB polarization so¬ 
lutions lead to consistent lensing 5-mode power spectrum mea¬ 
surements within 1 (T over the entire probed multipole range. In 
addition, the;^^ relative to the 2015 Planck ACDM base model 
(Planck Collaboration XIII 2015) of the lensing 5-mode band 
powers in eight multipole bins are 12.19, 8.99, 7.96, and 6.98, 
with corresponding probability-to-exceed values of 14, 34, 44, 
and 54 % for Coimnander, NILC, SEVEM, and SMICA respectively, 
which indicate the agreement of the Planck lensing B-mode 
signal with theoretical expectations. We measure amplitude fits 
with respect to the Planck base model of 


Aglens — 1.03 +0.11 
Allens = 0.97 ± 0.09 
Aglens = 1.02 ±0.10 
Allens = 0.97 + 0.08 


(Coimnander), 

(NILC), 

(SEVEM), 

(SMICA), 


corresponding to 10, 11, 10, and \2cr detections of lensing 
5 modes for Coimnander, NILC, SEVEM, and SMICA, respec¬ 
tively. 


9. Summary and recommendations 

We now summarize the results, and provide a critical analysis of 
the applicability of the derived maps for cosmological purposes. 

Starting with the temperature case, we have shown that the 
four CMB maps are in excellent agreement overall. The ampli¬ 
tudes of the pairwise difference maps are smaller than 5 yuK over 
most of the sky on large angular scales, and the high-^ power 
spectra agree to ~ 1 cr. Correspondingly, the differences with re¬ 
spect to the Planck 2013 maps are typically smaller than lOyuK 


over most of the sky, and the morphology of the differences is 
well understood in terms of improved treatment of systematic 
errors in the 2015 analysis. We conclude that the Planck 2015 
temperature maps provide a more accurate picture of the CMB 
sky than the 2013 temperature maps, having both lower noise 
and lower levels of systematics, and we expect these maps to find 
the same cosmological applications as the previous generation of 
maps. However, we emphasize that these maps are not cleaned 
of high-^ foregrounds, such as extragalactic point source or SZ 
emission. Residuals lead to biases in cosmological parameters at 
2 (T level, beyond I ^ 2000. Cosmological analyses using small 
angular scales must therefore take care to marginalize over such 
foregrounds as appropriate. 

For polarization, the situation is significantly more compli¬ 
cated due to two different problems. First, a low level of resid¬ 
ual systematics on large angular scales prevents a faithful CMB 
polarization reconstruction on multipoles i < 20. These modes 
are therefore removed by high-pass filtering in the current maps. 
Any cosmological analysis of these maps must take into ac¬ 
count the corresponding transfer function in order to avoid bi¬ 
ases. Second, due to the current noise mismatch between the 
FFP8 simulations and the data, we strongly caution against us¬ 
ing the polarization maps provided here for any cosmological 
analysis that depends sensitively on the assumed noise level. 
Nevertheless, the maps should prove useful for a number of other 
applications that do not require detailed noise simulations, for in¬ 
stance estimation of cross-spectra or cross-correlations, or stack¬ 
ing analyses, and we therefore release the maps to the public 
despite these limitations. 

Considering the four component separation methods in 
greater detail, the results may be distinguished according to two 
criteria, data selection and basis functions. While Coimnander 
performs data selection at the detector (or detector set) map 
level, rejecting potentially problematic maps, the other three 
methods employ frequency channel maps, and thereby maximize 
the signal-to-noise ratio. Fikewise, while Coimnander and SEVEM 
perform their analyses in pixel space, NILC and SMICA perform 
all operations in harmonic space. These distinctions can explain 
many of the qualitative differences discussed in the previous sec¬ 
tions. 

We make the following recommendations regarding the use 
of the four maps. First and foremost, we strongly recommend 
that any cosmological analysis based on these maps consider all 
four maps in parallel in order to assess the impact of specific 
choices of implementation and modelling. To be considered ro¬ 
bust, no results should depend strongly on the specifics of a given 
component separation algorithm. Considering specific details, 
we generally consider Coimnander to be the preferred solution 
on large and intermediate angular scales, due to its somewhat 
lower large-scale effective polarization noise (Table 1), weaker 
cross-correlation with the high-pass filtered WMAP K-Ka band 
synchrotron template (Table 2), lower A-point correlation func¬ 
tion fluctuations (Fig. 18 and Table 4), and weaker cosmologi¬ 
cal parameter dependence on ^max, suggesting less internal ten¬ 
sion between low, intermediate, and high multipoles (Fig. 15). 
In addition, the method is able to propagate uncertainties from 
the input maps to final products by means of Monte Carlo sam¬ 
ples drawn from the full posterior. For these reasons, we adopt 
the Coimnander solution for the low-^ Planck 2015 temperature 
likelihood (Planck Collaboration XI 2015). 

However, at high multipoles the Coimnander solution ex¬ 
hibits a significantly higher effective point source amplitude than 
the other three maps, due to the exclusion of frequencies be¬ 
low 217 GHz. For temperature, the lowest residual high-^ fore- 
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grounds are instead obtained by SMICA, as shown in Fig. 14. 
As a result, as in 2013, we confirm our preference for the SMICA 
map for analyses that require full-resolution observations in tem¬ 
perature, such as /jqL (fii'st three rows in Table 5 and Planck 
Collaboration XVII 2015) or lensing reconstruction (Fig. 19 
and Planck Collaboration XV 2015). SEVEM is also a very good 
choice for temperature, providing the map with the lowest level 
of noise at a wide range of scales as well as a smooth noise 
power spectrum, as measured by the HMHD maps (see Table 1 
and Fig. 14). It also performs equally well as SMICA with regard 
to the estimation of /^l- 

In polarization, NILC and SMICA perform equivalently at 
high multipoles (Fig. 14). The NILC polarization maps yield 
measurements of which are most consistent with zero 
(Table 5). The NILC and SMICA analyses also provide an ef¬ 
fective mapping of the weights of the Planck frequencies in the 
needlet/harmonic domains: a given weight tends to be high if a 
given frequency channel, in a given band in the harmonic do¬ 
main, is relevant for foreground cleaning; on the other hand, 
the higher the statistical noise at a given frequency, the lower 
the associated weight. For NILC and SMICA, the weights for T, 
E, and B modes are shown in Figs. B.2 and D.l, respectively. 
SEVEM provides the most stability with respect to the effect of 
bright point sources in polarization, (see Section 7.1). This is 
due to the inpainting procedure applied to these sources, which 
significantly reduces the effect of this contaminant. Therefore, 
SEVEM could be a suitable choice for those analyses in polariza¬ 
tion which cannot easily deal with the presence of point source 
holes in a mask. 

Finally, we note that the SEVEM approach is unique in its abil¬ 
ity to provide independent CMB estimates in a number of fre¬ 
quency channels. For analyses that benefit significantly from, or 
even require, such information, SEVEM is the only meaningful 
choice. Specific examples include various isotropy estimators 
(Planck Collaboration XVI 2015), the integrated Sachs-Wolfe 
stacking analysis (Planck Collaboration XXI 2015), relativis¬ 
tic boosting (Planck Collaboration XXVII 2014), and Rayleigh 
scattering analyses (Lewis 2013). 


10. Conclusions 

We have presented four different foreground-reduced CMB 
maps in both temperature and polarization derived from the 
Planck 2015 observations. These maps are based on the full 
Planck data, including a total of 50 months of LFI observa¬ 
tions and 29 months of HFI observations. The temperature com¬ 
ponent of these maps represents the most accurate description 
of the CMB intensity sky published to date. In the polariza¬ 
tion component, the characteristic £’-mode signal expected in a 
standard ACDM model is easily discernible over the full sky. 
Corresponding astrophysical foreground products are described 
in Planck Collaboration X (2015). 

The CMB maps presented here are the direct result of 
the detailed analyses of systematic errors described in Planck 
Collaboration VI (2015) and Planck Collaboration VIII (2015), 
which led to an effective reduction of systematic errors by almost 
two orders of magnitude in power on large angular scales in po¬ 
larization compared to 2013. However, despite these improve¬ 
ments, the polarization systematic errors in the Planck 2015 
data set are not yet negligible in several frequencies on the very 
largest scales. Multipoles below £ <20 are therefore suppressed 
by low-pass filtering in the current component-separated CMB 
maps. 


Additionally, as already noted in Planck Collaboration XII 
(2015), we observe a mismatch in the effective noise ampli¬ 
tude of 10-20 % when comparing the latest generation of Planck 
simulations (FFP8) with the data. Considering that the current 
temperature sky maps are signal-dominated up to ^ 2000, 

this noise mismatch is of little practical importance for cos¬ 
mological analyses based on temperature observations except 
on the very smallest scales. As in 2013, the temperature maps 
presented here are therefore used for a wide range of impor¬ 
tant applications, including large-scale temperature likelihood 
estimation (Planck Collaboration XI 2015), gravitational lens¬ 
ing (Planck Collaboration XV 2015), studies of isotropy and 
statistics (Planck Collaboration XVI 2015), primordial non- 
Gaussianity (Planck Collaboration XVII 2015), and non-trivial 
cosmological topologies (Planck Collaboration XVIII 2015). 
For high-^ power spectrum and likelihood estimation, we rec¬ 
ommend the cross-spectrum based methods described in Planck 
Collaboration XI (2015), primarily due to difficulties in estab¬ 
lishing sufficiently accurate models of unresolved extra-galactic 
high-^ foregrounds for the maps presented here. Cosmological 
parameters derived by temperature power spectra from these 
maps have been compared with the results of the Planck 2015 
likelihood analysis (Planck Collaboration XI 2015) and found to 
agree at the Icr level. 

For polarization, the noise mismatch is not negligible, and 
we therefore do not yet recommend using the provided maps for 
cosmological studies that require a highly accurate noise model. 
The polarization maps presented here may still be very useful 
for many important cosmological applications, including cross¬ 
correlation and cross-spectrum based analyses, and we therefore 
release the maps despite the current noise mismatch. Analysis 
of the higher order statistics of these maps has been performed 
within the current framework of precision assessment and pre¬ 
sented in this paper. The bispectrum analysis, including E mode 
polarization, applied to the results of the four component sepa¬ 
ration methods, gives evidence of a vanishing non-Gaussian sig¬ 
nal for three geometrical configurations, namely local, orthogo¬ 
nal, and equilateral. The 5-modes derived from the current po¬ 
larization maps have been cross-correlated with the predicted 
lensing 5-modes from the measured E signal and the lensing 
potential measured independently from temperature in Planck 
Collaboration XV (2015); the result is found to be in excellent 
agreement among the four component separation methods, as 
well as with the prediction of the Planck best-fit cosmology. 

On the basis of these encouraging results, intense work is on¬ 
going to reduce the large-scale polarization systematics to negli¬ 
gible levels, as well as to resolve the noise simulation mismatch, 
and good progress is being made. Updated products will be pub¬ 
lished as soon as this work has reached a successful completion. 
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Appendix A: Bayesian parametric fitting 

Coimnander (Eriksen et al. 2004, 2008) fits a physical model to a 
set of observations within a standard Bayesian parametric frame¬ 
work, defined by a set of explicit physical parameters and priors. 
The code can be operated in two modes, either employing Gibbs 
sampling to map out the full parameter posterior, or using iter¬ 
ative non-linear searches to derive the maximum-likelihood so¬ 
lution; the only implementational difference between the two is 
whether to sample from or maximize the conditional posterior in 
each Gibbs step. All maps presented in this paper are derived in 
the maximum-likelihood mode, and uncertainties are evaluated 
through simulations. 

Coimnander forms the core of the Planck 2015 foreground- 
targeted diffuse component separation efforts, and the cor¬ 
responding results are described in full detail in Planck 
Collaboration X (2015). In this section we only summarize the 
most relevant steps for CMB-oriented analysis. 



Fig.A.l. Commander processing masks for temperature {top) 
and polarization {bottom). For temperature, the different shades 
of grey correspond to different angular resolutions, ranging from 
5' (light grey) through 7.5' (dark grey) to 40' FWHM (black). 
For polarization, the same mask is used for both 10' and 40' 
FWHM resolution. 


A.1. Intensity 

In 2013, the Commander CMB temperature solution was de¬ 
rived using only the seven lowest Planck frequency maps be¬ 
tween 30 and 353 GHz, adopting a simple four-component sig¬ 
nal model, including CMB, CO, greybody thermal dust and a 
single power-law low-frequency component. The only instru¬ 
mental parameters included in the analysis were monopole and 
dipole corrections. In the current release, significantly more data 
are included in the analysis, and the astrophysical and instru¬ 
ment models have been expanded to account for more effects. 
Specifically, a total of 32 maps are considered in the analysis, in¬ 
cluding 21 Planck detector and detector set maps, 10 WMAP dif¬ 
ferencing assembly maps, and a 408 MHz low-frequency survey 
map (Haslam et al. 1982). This wide frequency range allows us 
to fit separately for synchrotron, free-free, spinning dust, thermal 
dust, and CMB, as well individual CO transitions at 115, 230 and 
353 GHz, a common line emission component in the 94 and 100 
GHz channels (primarily HCN), and thermal Sunyaev Zeldovich 
emission for the Coma and Virgo clusters. On the instrumental 
side, we now fit for both calibration and bandpass uncertainties, 
in addition to monopoles and dipoles (Planck Collaboration X 
2015). 

A.2. Polarization 

The Coimnander CMB polarization map is derived in an analo¬ 
gous manner to the temperature solution, but relying on Planck 
observations alone. At low resolution, we derive a 40' map from 
all frequencies between 30 and 353 GHz, including CMB, syn¬ 
chrotron and thermal dust in the signal model. At high resolu- 
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Fig.A.2. Multipole moment weights used for multi-resolution 
hybridization in the Commander CMB map, as described by 
Eq. A.l. 



tion, we derive a 10' map from frequencies between 100 and 
353 GHz, including only CMB and thermal dust. 

As discussed in Planck Collaboration X (2015) several spec¬ 
tral index models have been explored for polarized synchrotron 
and thermal dust emission, including 1) spatially constant, 2) 
smoothed over large angular scales, and 3) based on the temper¬ 
ature model. The main conclusion, however, is that the current 
polarization data have too low signal-to-noise ratio for both syn¬ 
chrotron and thermal dust indices to discriminate between the 
three models at a statistically significant level. Allowing such ad¬ 
ditional degrees of freedom only increases the degeneracies be¬ 
tween the various components without improving the overall fit, 
and also leaves the solution sensitive to large-scale residual sys- 
tematics. For now, we therefore adopt the temperature-derived 
spectral parameters for both synchrotron and thermal dust for 
the primary Commander CMB polarization map. 

A.3. Hybridization of multi-resolution sky maps 

As currently implemented through pixelized fits, the Commander 
analysis requires uniform angular resolution across frequencies 
in order to estimate spectral parameters correctly. This implies 
that all frequency maps must be smoothed to the resolution of 
the lowest resolution channel before analysis. 

In the 2013 Planck release, this problem was partially solved 
by first determining spectral parameters at low resolution us¬ 
ing Commander, as described above, and then solving for the 
component amplitudes from full resolution data using a so- 
called Ruler step, leading to the Commander-Ruler hybrid. 
In the current release, we adopt a Commander-only multi-stage 
approach, in which the system is solved at four different an¬ 
gular resolutions, using different subsets of the data, but each 
with internally coherent angular resolutions. Explicitly, we first 
solve for the full parameter set at 1° resolution with temper¬ 
ature data only, combining Planck observations with external 
data (WMAP and the 408 MHz Haslam observations; see Planck 
Collaboration X 2015 for full details). We then fix the global 
parameters (monopoles, dipoles, calibration and bandpass cor¬ 
rections), eliminate the external data, and solve again using 


Fig.A.3. 5° X 5° zoom-in of the multi-resolution contributions 
to the Commander hybrid CMB map from the 40' {top left), 
7.'5 {top right) and ^^5' {bottom left) solutions, centered on the 
South Galactic Pole. The hybrid map is shown in the bottom 
right panel. 


only Planck observations at 40' resolution, while at the same 
time simplifying the low-frequency foreground model. Next, we 
eliminate all frequencies below 143 GHz, and solve for CMB, 
CO and thermal dust at 7.'5 resolution, before finally eliminat¬ 
ing also the 143 GHz channel and the CO component, and solve 
only for CMB and thermal dust at the native resolution of the 
217 GHz channel, roughly 4.'8. 

Thus, a series of CMB estimates is established, ranging from 
low resolution derived within a complete foreground model, 
to high resolution derived within a greatly reduced foreground 
model. The basis for the low-^ likelihood is the lowest resolution 
solution, implementing the most complete foreground model and 
exploiting both Planck and external data. No further processing 
is required to cover angular scales up to multipoles oil < 250. 

For high-resolution temperature CMB analysis, we hybridize 
the three Planck-orAy solutions, ranging between 40 and 5', into 
a single map as follows. We first define a mask for the low- 
resolution level by simple thresholding. For the second level, 
we also threshold the corresponding;^^ map, but we additionally 
exclude point sources, and we require that all pixels excluded 
by the lower resolution maks are excluded by the higher resolu¬ 
tion mask. This is repeated for the third and highest level, but in 
addition we also exclude all pixels with a CO amplitude larger 
than 0.5 Kkm/s (see Planck Collaboration X 2015), as CO is no 
longer included in the foreground model. The resulting masks 
admits 98.4, 95.0 and 82.4% of the sky, respectively, and are 
shown in the top panel of Fig. A.l. For polarization, we con¬ 
struct a similar mask from the product of the thresholded low- 
resolution and CO maps, and the same mask is applied at 
both 40 and 10' FWHM. The resulting mask admits 83.6% of 
the sky, and is shown in the bottom panel of Fig. A.l. 

Next, at each level the corresponding masked regions are re¬ 
placed with a constrained Gaussian realization (Eriksen et al. 
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2004), as drawn from P{a^^\C£, d), before co-adding the three 
maps in harmonic space using the following cosine apodization 
scheme, 


^fullres 


^40 


,., 40/1 ,., 7 . 5 \ J.5 , ,., 7 . 5^5 

Wn (i-w. ) -\-w, ^ 




adopting cosine apodization weights given by 


(A.l) 


w 


.40 


W 


7.5 

b 


1 for ^<200 

\[1 - cos( 7 r 3 f^)] for 200 < ^ < 300 

1 for ^ < 1000 

\[1- cosCtty^^)] for 1000 < ^ < 1200 


(A.2) 

(A.3) 


as illustrated in Fig. A.2. Figure A.3 shows the individual con¬ 
tributions from the various resolutions, as well as the final sum. 

The result is a single full-resolution CMB map with variable 
effective sky fraction as a function of multipole. For analyses 
employing only multipoles below ^ < 200, a total of 98.4 % 
of the sky corresponds to direct CMB measurements, with the 
rest being filled with a proper Gaussian constrained realization. 
Therefore, for very low multipoles, for which the partially miss¬ 
ing modes can be nearly exactly reproduced by the high-latitude 
information, the full sky is in practice available for analysis. For 
slightly smaller scales, 98 % of the sky is available for direct 
analysis. Correspondingly, below € < 1000 a total of 95.0% 
of the sky represents direct measurements, while at higher mul¬ 
tipoles, only 82% of the sky corresponds to direct measure¬ 
ments. When using the 2015 Coininander CMB map for cosmo¬ 
logical analysis, it is therefore important to consider which an¬ 
gular scales are relevant before choosing the appropriate mask. 
However, using the most conservative mask is always safe at any 
angular scale. 

Analogous processing is performed for the Commander po¬ 
larization map, but employing only two resolution levels, namely 
40' and 10' FWHM. Cosine hybridization is performed between 
^ = 200 and 300, similar to the lowest two levels in the temper¬ 
ature analysis. 


Appendix B: Internal linear combination in needlet 
space 

The goal of NILC is to estimate the CMB from multifre¬ 
quency observations while minimizing the contamination from 
Galactic and extragalactic foregrounds, and instrumental noise. 
The method makes a linear combination of the data from the 
input maps with minimum variance on a frame of spherical 
wavelets called needlets (Narcowich et al. 2006). Due to their 
unique properties, needlets allow localized filtering in both pixel 
space and harmonic space. Localization in pixel space allows the 
weights of the linear combination to adapt to local conditions 
of foreground contamination and noise, whereas localization in 
harmonic space allows the method to favour foreground rejec¬ 
tion on large scales and noise rejection on small scales. Needlets 
permit the weights to vary smoothly on large scales and rapidly 
on small scales, which is not possible by cutting the sky in zones 
prior to processing (Delabrouille et al. 2009b). 

The NILC pipeline (Basak & Delabrouille 2012, 2013) is ap¬ 
plicable to scalar fields on the sphere, hence we work separately 
on maps of temperature and the E and B modes of polarization. 
The decomposition of input polarization maps into E and B is 
done on the full sky. At the end, the CMB Q and U maps are 
reconstructed from the E and B maps. 


Prior to applying NILC, all of the input maps are convolved 
or deconvolved in harmonic space to a common resolution cor¬ 
responding to a Gaussian beam of 5' FWHM. Each map is then 
decomposed into a set of needlet coefficients. For each scale 7 , 
needlet coefficients of a given map are stored in the form of a sin¬ 
gle HEALPix map. The filters used to compute filtered maps 
are shaped as follows 
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For each scale 7 , the filter has compact support between the mul¬ 
tipoles /A^ and /^ax with a peak at /^^^ (see table B.l and figure 
B.l). The needlet coefficients are computed from these filtered 
maps on HEALPix pixels with Aside equal to the smallest power 
of 2 larger than /max/ 2 . 


Table B.l. List of needlet bands used in the NILC analysis. 


Band 


^peak 

/ 

^max 

Aside 

;• = 1 . 

0 

0 

100 

64 

2 . 

0 

100 

200 

128 

3 . 

100 

200 

300 

256 

4 . 

200 

300 

400 

256 

5 . 

300 

400 

600 

512 

6 . 

400 

600 

800 

512 

7 . 

600 

800 

1000 

512 

8 . 

800 

1000 

1400 

1024 

9 . 

1000 

1400 

1800 

1024 

10. 

1400 

1800 

2200 

2048 

11 . 

1800 

2200 

2800 

2048 

12. 

2200 

2800 

3400 

2048 

13. 

2800 

3400 

4000 

2048 


In order to show the contribution of the various frequency 
channels to the final CMB map at different needlet bands, we 
compute the full sky average of needlet weights for each fre¬ 
quency channel and needlet band. Figure B.2 shows that the 
most of contribution to the reconstructed CMB maps comes 
from the 143 GHz and 217 GHz channels. In the low i needlet 
bands, the contribution from 143 GHz is large compared to that 
from 217 GHz. However, due to better angular resolution, the 
217 GHz channel contributes more than the 143 GHz channel in 
the highest I needlet bands. In the intermediate needlet bands, 
the contributions from these two channels are comparable. It is 
interesting to note the contribution from the LFI in the lowest I 
bands. In intensity, the 70 GHz channel serves mostly for fore¬ 
ground removal, while for polarization it contributes positively 
to the CMB solution very similarly between E and B. We stress 
again that the lowest i modes have been filtered out in the re¬ 
sults presented here, and therefore the low I results will require 
further investigation. 
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Fig.B.l. Needlet bands used in the analysis. The solid black line 
shows the normalization of the needlet bands, that is, the total 
filter applied to the original map after needlet decomposition and 
synthesis of the output map from needlet coefficients. 


B.1. Masking 

The confidence masks for NILC for intensity and polarization 
have been generated following a procedure similar to that used 
by SMICA, but adopting a different parameterization. 

For intensity, the NILC CMB map is filtered through a spec¬ 
tral window 


/(/) = exp 


--1 

/- i7ooy 

_ A 

200 ) 


(B.l) 


The result is then squared and smoothed with a Gaussian circular 
beam with FWHM 120'. The variance map obtained in this way 
is then corrected for the noise contribution by subtracting the 
variance map for the noise obtained in the same way from the 
NILC HRHD map. The confidence map is obtained by thresh¬ 
olding the noise-corrected variance map at 73.5 juK^. For polar¬ 
ization, the polarized intensity P = ^JQ^ + is obtained from 
the NILC outputs, and is filtered through a circularly symmet¬ 
ric Gaussian window function of FWHM 30'. The result is then 
squared and smoothed with a Gaussian of FWHM 210'. The 
resulting variance map is corrected for the noise contribution 
by following the same procedure used for intensity. The confi¬ 
dence map is obtained by thresholding at 6.75 yuK^. The resulting 
masks are shown in Fig. B.3. 

Appendix C: Template fitting 

The SEVEN method (Leach et al. 2008; Fernandez-Cobos et al. 
2012) aims to produce clean CMB maps for several frequency 
channels by using internal template fitting. The templates are 
constructed from the Planck data, typically as the subtraction of 
two close Planck frequency channels to remove the CMB signal. 
The resolution of the maps are equalized before subtraction. The 
cleaning is achieved simply by subtracting a linear combination 
of the templates tj{x) from the data, with coefficients aj obtained 
by minimizing the variance outside a given mask. 


nt 

Tc(x, v) = d(x, v) - y] 
j=i 


(C.l) 



Frequency (in GHz) 



Frequency (in GHz) 

Fig.B.2. Full-sky average of needlet weights for different fre¬ 
quency channels and needlet bands. From top to bottom, the 
panels show results for temperature, E, and B modes. 


where nt is the number of templates used, and Tdx, v) and d(x, v) 
correspond to the cleaned and raw maps at frequency v, respec¬ 
tively. The same expression applies for T, 2, or U. 

The cleaned frequency maps are then combined in harmonic 
space, taking into account the noise level, resolution, and, for 
temperature, a rough estimate of the foreground residuals of each 
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NILC T 




Fig.B.3. NILC masks for temperature {top) and polarization 
{bottom). 


three templates, the first channel in the subtraction is smoothed 
with the beam of the second map and vice versa. For the 
857 GHz template, we simply filter the map with the 545 GHz 
beam (this is done for comparison with the 857-545 template 
from the 2013 pipeline, where the 857 GHz map was smoothed 
by the 545 GHz beam.) Note that the coefficients used to multi¬ 
ply this template are typically ~ 10“"^, so the level of CMB signal 
introduced by this template in the final cleaned map is negligi¬ 
ble. We take advantage of this to drop subtraction of the 545 GHz 
map, as was done in 2013 nominally to remove the CMB signal. 
This simplifies the method and also reduces the noise in this tem¬ 
plate. 

The 100, 143, and 217 GHz maps are then cleaned by sub¬ 
tracting a linear combination of the four templates. The coef¬ 
ficients of the linear combination (Table C.l) are obtained by 
minimizing the variance outside an analysis mask. The main dif¬ 
ference with respect to the 2013 release is that we have used 
the same coefficients for the whole sky (instead of dividing it in 
two regions), since this simplifies the procedure without affect¬ 
ing the quality of the reconstruction (other than on those pix¬ 
els very close to the Galactic centre that need to be masked in 
any case). This analysis mask covers the 1 % of the sky with the 
brightest emission, as well as sources detected at all frequency 
channels. Once the maps are cleaned, each is inpainted on the 
source positions detected at that (raw) channel. Then, the MHW 
algorithm is run again, now on the cleaned maps. A relatively 
small number of new sources are found and are also inpainted at 
each channel. The resolution of the cleaned map is the same as 
that of the raw map. Note that no assumptions about the noise or 
foregrounds are made in order to construct the single-frequency 
cleaned maps. 


cleaned channel, to produce a final CMB map at the required 
resolution. 

C.1. Implementation for temperature 

For temperature, we followed a similar procedure to that for 
the Planck 2013 release: the 100, 143, and 217 GHz maps are 
cleaned using four templates constructed from the six remaining 
frequency channels. A few differences have been implemented 
in the current pipeline with respect to the previous work: the use 
of a single coefficient over the whole sky for each template (in¬ 
stead of defining two regions); the use of inpainting to reduce 
contamination from sources; the use of the 857 GHz channel as 
a template (instead of 857-545). Note that the other three tem¬ 
plates (30-44, 44-70, and 353-217) are the same as for the pre¬ 
vious release. 

The six frequency channels used to construct templates (30 
to 70 GHz and 353 to 857 GHz) are inpainted at the position of 
sources detected by the Mexican hat wavelet (MHW) algorithm 
(Planck Collaboration XXVI 2015). The size of the holes to be 
inpainted is determined by taking into account the beam size of 
the channel as well as the flux density of each source. We do 
a simple diffuse inpainting, which fills one pixel with the mean 
value of the neighbouring pixels in an iterative way. To avoid 
inconsistencies when subtracting two channels, each map is in¬ 
painted on the sources detected in both that map and on the other 
map used to construct the template. For example, to construct the 
(30-44) template, both maps are inpainted in the positions of the 
sources detected at 30 and 44 GHz. This reduces significantly the 
contamination from compact sources in the templates. 

Once they have been inpainted, the maps are smoothed to 
a common resolution and then subtracted. To construct the first 


Table C.l. Linear coefficients, aj, of the templates used to clean 
individual frequency maps with SEVEN for temperature. 


Coefficients aj 


Template 100 GHz 143 GHz 217 GHz 


30-44 . -6.38 X 10“^ 2.84 x 10"^ -1.41 x 10-^ 

44-70 . 3.53 x 10-1 1.33 x lO"! 3.82 x lO"! 

545-353 . 4.34 x 10-^ 6.56 x 10-^ 1.79 x 10"^ 

857 . -3.63 x 10-5 -5.66x10-5 -1.18 xlO-^ 


Finally, the SEVEN CMB map is constructed by combining 
the cleaned and inpainted 143 and 217 GHz maps. In the com¬ 
bination, the maps are weighted taking into account the noise, 
resolution, and a rough estimation (obtained from realistic simu¬ 
lations) of the foreground residuals at each map. The resolution 
of this map corresponds to a Gaussian beam of FWHM 5' and 
HEALPix resolution Aside = 2048 with a maximum multipole 
^max = 4000. 

The same procedure, including the full inpainting process 
of point sources, is applied when running the pipeline on FFP8 
simulations. 

C.2. Implementation for polarization 

To clean the polarization maps, a procedure similar to the one 
used for the temperature maps is applied to the Q and U maps 
independently. However, given that narrower frequency cover¬ 
age is available for polarization, the templates and maps to be 
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cleaned are different. In particular, we clean the 70, 100, and 
143 GHz maps using three templates. Since the signal-to-noise 
ratio is lower, at 100 and 143 GHz the clean maps are produced 
at A^side = 1024, with resolution corresponding to a Gaussian 
beam of FWHM 10' and a maximum multipole Imax = 3071. At 
70 GHz, the map is produced at its native resolution. 

The first step of the pipeline is to inpaint the positions of the 
sources detected using the MHW algorithm in those channels 
that will be used to construct templates. As for the temperature 
case, for a given template constructed as the difference of two 
frequency channels, the inpainting is performed for all of the 
sources detected in any of the channels. Note that the inpainting 
is performed in the frequency maps at their native resolution. 

These inpainted maps are then used to construct a total of 
four templates. To trace the synchrotron emission, we construct 
a 30-44 template. For the dust emission, the following templates 
are used: 353-217 (smoothed to 10' resolution), 217-143 (used 
to clean 70 and 100 GHz) and 217-100 (to clean 143 GHz). 
These last two templates are constructed at 1° resolution, since 
an additional smoothing becomes necessary in order to increase 
the signal-to-noise ratio of the template. Otherwise, the esti¬ 
mated coefficients are driven by the noise and the cleaning is less 
efficient. Since fewer frequency channels are available in polar¬ 
ization, it becomes necessary to use the maps to be cleaned as 
part of one of the templates. Therefore the 100 GHz map is used 
to clean the 143 GHz frequency channel and vice versa, making 
the clean maps less independent than in the temperature case. 

These templates are then used to clean the non-inpainted 70 
(at its native resolution), 100 (at 10' resolution), and 143 GHz 
maps (also at 10'). The corresponding linear coefficients (listed 
in Tables C.2 and C.3) are estimated independently for Q and 
U by minimizing the variance of the clean maps outside a mask 
that covers compact sources and the 3 % of sky with the brightest 
Galactic emission. Once the maps have been cleaned, inpainting 
of the sources detected at each map is carried out. The size of the 
holes to be inpainted takes into account the additional smooth¬ 
ing of the 100 and 143 GHz maps. As for temperature, the same 
inpainting processing is applied to the point sources positions 
when running the pipeline on the FFP8 simulations. The 100 
and 143 GHz clean maps are then combined in harmonic space, 
using a full-sky E and B decomposition, to produce the final 
CMB maps for the Q and U components with a Gaussian beam 
of FWHM 10' and HEALPix resolution Aside = 1024. Each map 
is weighted taking into account its noise level at each multipole. 


Table C.2. Linear coefficients, aj, of the templates used to clean 
individual frequency maps with SEVEN for Q. 


Coefficients aj 


Template 70 GHz 100 GHz 143 GHz 


30-44 . 2.66x 10-2 1.00 xlO-^ 5.72x10-2 

217-143. 8.01x10-2 1.24x10-1 

217-100 . ... ... 1.78x 10-1 

353-217 . 4.26 x 10-^ 1.14 x IO -2 2.50 x IO -2 


Before applying the post-processing high-pass filter to the 
cleaned Q and U polarization maps, we inpaint the region with 
the brightest Galactic residuals (5 % of the sky) with the same 
simple algorithm used for point source holes. This is done to 


Table C.3. Linear coefficients, aj, of the templates used to clean 
individual frequency maps with SEVEN for U. 


Coefficients aj 


Template 70 GHz 100 GHz 143 GHz 


30-44 . 2.98 x 10-2 1.47x 10-2 7.97 x IO -2 

217-143 . 4.08 x 10-2 1.53 x lO-i 

217-100 . ... ... 1.66x 10-1 

353-217 . 6.76 x IO -2 0.66 x IO -2 2.41 x IO -2 


avoid introducing ringing around the Galactic centre when the 
maps are filtered. 

In addition, E and B modes maps are constructed from the 
clean Q and U maps. In this case, prior to performing the decom¬ 
position, the region of the Q and U maps defined by the SEVEN 
confidence mask is filled with a Gaussian-constrained realization 
(Eriksen et al. 2004). 

C.3. Masks 

In temperature, the SEVEN confidence mask is produced by 
thresholding differences obtained between different CMB re¬ 
constructions. In particular, we construct the difference map be¬ 
tween the clean 217 and 143 GHz maps at a resolution of FWHM 
30' and Aside = 256. The brightest pixels of this map (and 
their direct neighbours) are successively removed from the CMB 
combined map (at the same resolution) and the dispersion of 
the CMB combined map calculated. If a sufficient number of 
pixels is removed, the dispersion of the CMB map goes down 
and reaches a plateau, indicating that convergence has been 
achieved. The removed pixels constitute the mask. This mask 
is then smoothed with a Gaussian beam of 1° to avoid sharp 
edges, and upgraded to full resolution. The same procedure is 
repeated for other two difference maps: the clean 143-100 map 
and the difference of two clean CMB combined maps, whose 
linear coefficients have been obtained by minimizing the vari¬ 
ance outside two different masks. Finally, the three masks are 
multiplied in order to produce our final confidence mask, which 
leaves a suitable sky fraction of approximately 85 %. Residual 
monopoles and dipoles outside this mask are subtracted for the 
single frequency and combined cleaned maps. 

For polarization, the clean combined map is downgraded to 
a resolution of FWHM 90' and Aside = 128. The dispersion at 
each pixel is estimated from a circle centered in the considered 
pixel. Those pixels with a dispersion above a given threshold are 
included in the mask, which is then smoothed with a Gaussian 
beam of FWHM 90' to avoid sharp edges and upgraded to the 
required Aside- Finally, this mask is multiplied by a mask cus¬ 
tomized to cover the CO emission, in order to discard those pix¬ 
els contaminated by this foreground component due to the band¬ 
pass leakage. An additional 1 % of the sky is added to the final 
mask to remove those pixels most affected by the high-pass fil¬ 
tering that is subsequently applied to the cleaned maps. The fi¬ 
nal mask allows for a useful fraction of the sky of approximately 
80 %. The 2015 SEVEN masks are shown in Fig. C.l. 

Appendix D: Spectral matching 

Spectral Matching Independent Component Analysis (SNICA) is 
a semi-blind component separation algorithm which operates in 
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SEVEM T 




Fig. C.l. SEVEM masks in temperature {top) and polarization 
{bottom). 


harmonic space. CMB maps are synthesized from spherical har¬ 
monic coefficients, obtained as weighted linear combina¬ 
tions of the coefficients of A/chan input maps, 




(D.l) 


where is the Auchan x 1 vector of the spherical harmonic co¬ 
efficients of the input maps and is the A/chan x 1 vector of 
weights. The spectral weights used to produce the SMICA maps 
are designed to minimize the total foreground and noise contam¬ 
ination at each multipole, under the constraint that the resulting 
map has a well defined effective beam window function, that of 
a Gaussian beam with 5' FWHM. 

In theory, such weights are given by 




a^R“^a ’ 


(D.2) 


where the Auchan x 1 vector a is the frequency spectrum of the 
CMB and the A/chan x Auchan spectral covariance matrix con¬ 
tains in its (/, j) entry the cross-power spectrum at multipole i be¬ 
tween the input frequency maps i and j, at 5' resolution. In prac¬ 
tice, the spectral covariance matrices must be estimated from the 
data. The natural sample estimate 




tm 


(D.3) 


can be used, possibly after some binning, to replace R^ in the 
weight formula (D.2). This works well at high ^ because a large 
number of modes are averaged in (D.3), so this estimate has 


low variance. At larger scales, however, it is necessary to con¬ 
strain the spectral covariance matrices in order to get reliable es¬ 
timates. For that purpose, SMICA uses a semi-parametric model 
of these matrices. 

The CMB, the foregrounds, and the noise are independent 
processes, so the spectral covariance matrices, after beam cor¬ 
rection, can be decomposed into 

Re = + Rf'’ + Rf(D.4) 

We assume that the noise is uncorrelated between frequency 
channels, therefore the noise contribution is diagonal. The 
signal parts are modelled as 

R^'"'’ = aa^Q, Rf‘^ = FP^F^ (D.5) 

where Q is the CMB power spectrum, matrix F is A/chan x d, and 
is a J X J positive definite matrix. A SMICA fit consists of fit¬ 
ting the model (D.4) and (D.5) to empirical covariance matrices 
R^ by minimizing the spectral matching criterion 

^(2^ -f 1) [trace(R^ + IndetR^j. 

The fitted values of R^ can be seen as regularized versions of 
their empirical counterparts R^, to be used in equation (D.2). If 
no constraints are imposed on matrices F and P^, except that the 
latter is positive definite, the foreground model is equivalent to 
d sky templates with arbitrary frequency spectra (represented by 
the columns of F) and arbitrary power spectra and correlations 
between them (represented by P^). Ultimately, the foreground 
contribution is controlled by a single parameter, the rank d of 
the foreground model. 

More details regarding the principles of SMICA can be found 
in Cardoso et al. (2008). We now describe its extension to polar¬ 
ization data. 

There are several options to extend SMICA to polarization. 
We choose to obtain the CMB Q and U maps through a joint 
processing of the E and B modes. The CMB U and Q maps are 
synthesized from their E and B modes, and obtained as 
linear combinations 



where the A/chan x 1 vectors and x^^ are the spherical har¬ 
monic coefficients of the input maps and is a 2Achan x 2 ma¬ 
trix of weights. The optimal weights are obtained as a simple 
generalization of the single-field case. 


W, 


(A^R^'a) ‘a^R^‘ with A = 


a 0 
0 a 


(D.7) 


where now is a 2Achan x 2Achan covariance matrix. The 
weights defined by (D.7) can be safely obtained at high £ by us¬ 
ing the sample covariance matrices, possibly after binning them. 
At low multipoles, some regularization via modelling is in or¬ 
der, as for the temperature analysis. We use a natural extension 
of temperature (D.5), in the form 


R^”'' = A 


f^EE f^EB 
(jEB (jEE 


A^ 


R 


fgd 


Fe 0 

P. 

F£ 0 

O 

to 


o 

to 


(D.8) 

(D.9) 
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where and are Auchan x d matrices and is a 2J x 2d 
matrix. For this release, no constraints are imposed on those ma¬ 
trices, except that each must be positive definite. Hence, as in 
temperature, SMICA fits a foreground model representing d po¬ 
larized templates with arbitrary frequency spectra, power spec¬ 
tra, and correlation. 

D.1. Implementation for temperature 

The production of the CMB temperature map is mostly un¬ 
changed from 2013 (see Planck Collaboration XII 2014 for de¬ 
tails). Here, we recall the 3-step fitting strategy adopted in 2013 
and mention the changes made for this release. 

- First fit: recalibration. A preliminary and independent SMICA 
fit is done with power spectra estimated over a clean fraction 
of the sky, including a as a free parameter. This can be un¬ 
derstood as a recalibration procedure; the estimated value of 
a is kept fixed in the following steps. The number of fore¬ 
ground templates, J, is now 5, whereas it was 4 in 2013. The 
30 GHz channel is not recalibrated, unlike in 2013. 

- Second fit: foreground emission. The foreground emission, 
matrix F, is estimated in a second SMICA fit with a kept fixed 
at the value found in the first step. Spectra are estimated over 
a large fraction of the sky (/sky = 97 %) and the fit is made 
over the range 4 < ^ < 150. For the data, this step is the same 
as in 2013. For FFP8, foreground emission is captured using 
7 templates, whereas 6 templates were used for the FFP6 
simulations in 2013. 

- Third fit: power spectra. The frequency spectra captured by 
vector a and matrix F are kept fixed at the values found in 
steps 1 and 2; we fit only the signal power spectra and 

and the frequency channel noise power spectra. Spectral 
covariance matrices are computed over 97 % of the sky and 
the fit includes all multipoles up to ^ = 1500. This step is the 
same as in 2013. 

The weights determined for the temperature maps are shown 
at the top of Fig. D.l. It should be noted that the CMB map 
is synthesized from spherical harmonic coefficients which have 
been set explicitly to zero for ^ = 0 and £ = 1. Therefore, the 
SMICA CMB map has no monopole or dipole components over 
the full sky. 

D.2. Implementation for polarization 

All Planck polarized channels are used to produce the polarized 
CMB maps. The process is easier in polarization than it is in 
temperature because less precision is required due to the lower 
signal-to-noise ratio and also because the foregrounds appear to 
have a simpler structure. In particular, we do not preprocess the 
frequency maps by subtracting or masking point sources as is 
done to the temperature maps. 

The SMICA fit in polarization is conducted with the same pa¬ 
rameters as in temperature, but with two differences. First, the 
recalibration step is omitted: we use the CMB frequency spec¬ 
trum (vector a) determined from temperature maps. Second, the 
foreground model comprises 6 polarized templates (as for tem¬ 
perature) but the matrices F e and F b are fitted in the second step 
over the multipole range 4 < ^ < 50. The weights determined 
for the polarization data are shown at the bottom of Fig. D.l. 

In order to mitigate spectral leakage from E io B and from 
the Galactic plane onto other regions of the sky, we do not com¬ 
pute spherical transforms directly from masked Q and U maps. 




Fig. D.l. SMICA weights for temperature {top) and polarization 
{bottom). For readibility, the values are shown for input maps 
in units of antenna temperature. The plot goes up to 7” ~ 3600, 
but the output maps are synthesized uses all multipoles up to 
7 = 4000. For polarization, the thick solid lines show the con¬ 
tribution of input E modes to the CMB E modes and the thick 
dashed lines show the same for the B modes. The thin lines, all 
close to zero, show “cross-contributions” of input E modes to 
the CMB B modes and vice versa. 

Instead, for each pair of input Q and U maps, we first produce 
full-sky E and B maps to which an apodized Galactic mask is 
applied. It is from those masked E and B maps that spherical har¬ 
monic coefficients are computed for the estimation of the spectra 
and cross-spectra going into and for the synthesis of the final 
CMB map. 

D.3. Masks 

Confidence masks, shown in Fig. D.2, are built using the pro¬ 
cedure described in Planck Collaboration XII (2014), with the 
following changes. For temperature, we apply a bandpass filter 
to the CMB map that is then squared and smoothed at 3?5 (it was 
2° in 2013). The confidence mask is obtained by thresholding the 
resulting map of local power. The threshold is determined by vi- 
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sual inspection to be 50yuK^ (it was 70/iK^ in 2013). The result¬ 
ing mask is then enlarged by multiplication by a Galactic mask 
covering 10 % of the sky. For polarization, the confidence mask 
is obtained by a procedure similar to temperature, but the CMB 
Q and U maps are low-pass filtered (rather than bandpassed) us¬ 
ing a Gaussian beam with 30' FWHM. They are then squared 
and smoothed to 3?5 resolution. Any area where the resulting 
P = a/G^ + map is above 5 is excluded from the confi¬ 
dence mask. In addition to that, we exclude pixels which are less 
than 7° away from the Galactic equator. 


SMICA T 



Fig.D.2. SMICA masks in temperature {top) and polarization 
{bottom). 


Appendix E: FFP8 Simulations 

In this appendix we provide a compendium of analyses eval¬ 
uated from the FFP8 data set, corresponding directly to those 
performed on the Planck 2015 data in the main text. 

E.1. CMB map differences 

Figures 4, 9 and 10 show pairwise difference maps between any 
of the Commander, NILC, SEVEM, and SMICA maps derived from 
the Planck 2015 data for each of the three Stokes parameters, /, 
G, and U. In Figs. E.1-E.3 we show the same evaluated from the 
FFP8 simulation set. 

Overall, the relative differences are of similar magnitude in 
the simulated data set as in the data, although they have some¬ 
what different dominant morphology. The high-latitude differ¬ 
ences are in general somewhat weaker in FFP8 than in the data, 
while the low-latitude differences are somewhat stronger. This 
is primarily due to the rather complicated foreground model 


adopted for thermal dust emission in the FFP8 simulations (see 
Planck Collaboration XII 2015 for details). In short, thermal 
dust emission is modelled in the FFP8 simulations by a sum 
of two greybody components, one of which has a temperature- 
dependent spectral index, P{T^. This was motivated by the re¬ 
sults presented by Planck Collaboration Int. XIV (2014); how¬ 
ever, as shown by the updated analysis in Planck Collaboration 
X (2015), there is no evidence for steepening from the Planck 
data alone. Only when the IRAS 100 pm data are included in the 
fit is any such effect seen. At such high frequencies the thermal 
dust physics is far more complicated, and the overall calibration 
problem more difficult. For the present discussion, it is sufficient 
to note that the thermal dust model adopted for the FFP8 simu¬ 
lations is significantly more complicated in terms of frequency 
dependence than the observed sky, and includes significant spec¬ 
tral curvature below 353 GHz that is not seen in the actual sky. 
None of the component separation codes accounts for this ad¬ 
ditional curvature, and this results in the strong features near 
the Galactic plane seen in some of the pairwise difference maps 
shown in Fig. 4. 

E.2. 1-point statistics of the total intensity FFP8 maps 

Next, we consider the 1-point statistics of the FFP8 fiducial map 
in temperature, using the FFP8-UT74 mask (Table E.l). These 
statistics were defined in Sect. 7.1, and corresponding polariza¬ 
tion results were shown in Table 3. We see that all the meth¬ 
ods are in good agreement with the Monte Carlo simulations at 
all resolutions considered, and the differences among codes are 
small. 

Table E.l. Lower tail probability in percent for the variance, 
skewness, and kurtosis for the FFP8 total intensity analysis at 
three different resolutions. The results have been obtained with 
the FFP8-UT74 mask. 



Lower-tail Probability 

Map 

Variance Skewness Kurtosis 


A^side - 2048 


Commander . 

81.1 

52.8 

37.0 

NILC. 

85.8 

47.8 

16.2 

SEVEM . 

76.1 

44.7 

30.6 

SMICA . 

82.4 

50.4 

38.2 

Nsiie = 256 

Commander . 

79.3 

57.5 

16.9 

NILC. 

85.0 

54.7 

12.9 

SEVEM . 

76.8 

55.5 

21.3 

SMICA . 

80.8 

55.2 

25.2 

Aside “ 64 

Commander . 

81.3 

74.0 

29.1 

NILC. 

85.0 

71.1 

23.8 

SEVEM . 

81.1 

72.5 

29.0 

SMICA . 

83.2 

71.1 

22.8 


E.3. The real-space N-point correlation functions for FFP8 
maps 

Finally, we present the real-space A-point correlation functions 
derived from the FFP8 simulations, complementing the analy- 
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Commander - NILC Commander - SEVEM 



-7.5 jiiK 7.5 


Fig.E.l. Pairwise difference maps between CMB temperature maps obtained on FFP8 simulations. Prior to differencing, the maps 
have been smoothed to 80 arcminutes FWHM and downgraded to NsUe = 128. 


sis shown in Sect. 7.2. As for the data, we analyse the high- 
pass filtered FFP8 CMB estimates downgraded to Aside = 64 and 
smoothed with a Gaussian kernel of 160' FWHM. For the polar¬ 
ization analysis we employ a low resolution version of the com¬ 
mon polarization mask, and for the temperature analysis we em¬ 
ploy a low resolution version of the common temperature mask. 

The A-point functions for the FFP8 simulations are shown in 
Figs. E.4 and E.5, plotted in terms of differences between the A- 
point functions for the high-pass filtered fiducial FFP8 map and 
the mean value estimated from 1000 Monte Carlo simulations. 
The probabilities of obtaining larger values for the statistic, 
compared to the Planck fiducial ACDM model, are tabulated in 
Table E.2. 

Overall, there is good agreement between the recovered 
FFP8 CMB estimates and the corresponding Monte Carlo en¬ 


sembles. The main outlier is NILC, for which the fluctuations 
are larger than expected for all three-point polarization A-point 
functions. 


E.4. CMB Power spectra and residuals from individual 
components 

Here we discuss the recovery of the power spectrum in the FFP8 
simulations, and assess the impact of foregrounds in the CMB 
solutions presented in this paper. CMB angular power spectra 
from the foreground cleaned CMB maps from the FFP8 are 
shown in Fig. E.6 compared to the input for TT, EE. The perfor¬ 
mance on simulations is comparable to that on the data, as may 
be seen by comparison with Fig. 14. 
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Commander - NILC 



Commander - SMICA 



NILC - SMICA 



Commander - SEVEM 



-1 


//K 


1 


Fig. E.2. Pairwise difference maps between Q maps obtained on FFP8 siulations. Smoothing and degrading as in Fig. 4. 


Figures E.7 and E.8 show the residual effect of combined 
individual components on the foreground-cleaned CMB maps. 
These plots give an estimate of how the residuals of the fore¬ 
ground components in the CMB solutions presented in this paper 
compare in power and are distributed in i. The FFP8 simulated 
skies of the individual components are propagated through the 
various component separation pipelines, and their power spec¬ 
tra are calculated using the FFP8-UT74 sky mask. The individ¬ 
ual components are those simulated in the FFP8 (Sec. 2.2 and 
Planck Collaboration XII (2015)) and include in TT dust, the 
far-infrared background (FIRB), unresolved sources (PS), and 
all the remaining components summed together (“other”). In 
EE, they include synchrotron, dust, noise, and remaining com¬ 
ponents. With the exception of the lowest multipoles, the domi¬ 
nant contaminant to the TT signal is noise. For EE, the residuals 


from dust dominate over the other foregrounds, but remain sub¬ 
dominant with respect to instrumental noise at all multipoles. 
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Commander - NILC 


Commander - SEVEM 



Fig,E.3. Pairwise difference maps between U maps obtained on FFP8 siulations. Smoothing and degrading as in Figs. 4,9. 





Fig.E.4. The difference between the A^-point functions and the corresponding means estimated from 1000 MC simulations. From 
left to right, results for the 2-point, pseudo-collapsed 3-point, equilateral 3-point and connected rhombic 4-point functions for the 
A^side = 64 FFP8 CMB temperature estimates. The black solid, red dot dot dot-dashed, orange dashed, green dot-dashed, and blue 
long dashed lines correspond to the true. Commander, NILC, SEVEM, and SMICA maps, respectively. The true CMB map was analysed 
with added noise corresponding to the SMICA component separation method. The shaded dark and light grey regions indicate the 
68 % and 95% confidence regions, respectively, estimated using SMICA simulations. See Sect. 7.2 for the definition of the separation 
angle 6. 
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Fig.E.5. The difference between the A^-point functions for the high-pass filtered A^side = 64 FFP8 CMB estimates and the corre¬ 
sponding means estimated from 1000 MC simulations. The Stokes parameters Qr and Ur were locally rotated so that the correlation 
functions are independent of coordinate frame. The first row shows results for the 2-point function, from left to right, TQr, QrQr, 
and QrUr. The second row shows results for the pseudo-collapsed 3-point function, from left to right, TTQr, TQrQr, QrQrUr, 
and UrUrUr, and the third row shows results for the equilateral 3-point function, from left to right, TTQr, TQrQr, QrQrUr, and 
UrUrUr. The black solid, red dot dot dot-dashed, orange dashed, green dot-dashed, and blue long dashed lines correspond to the 
true. Commander, NILC, SEVEN, and SMICA maps, respectively. The true CMB map was analysed with added noise corresponding 
to the SMICA component separation method. The shaded dark and light grey regions indicate the 68% and 95% confidence regions, 
respectively, estimated using SMICA simulations. See Sect. 7.2 for the definition of the separation angle 6. 
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Table E.2. Probability-to-exceed (PTE) in percent for the Ap¬ 
point correlation function statistic applied to the FFP8 sim¬ 
ulation at A^side = 64 for each of the four CMB codes, as shown 
in Figs. E.4 and E.5. For reference, the second column lists the 
corresponding probabilities for the true input map with added 
noise corresponding to the SMICA map. 



PTE [%] 

Correlation function 

Input Comm. NILC SEVEM SMICA 


A/^-point; see Fig. E.4 


2-pt . 

30.4 

33.4 

40.4 

33.3 

30.7 

Pseudo-coil. 3-pt .... 

11.1 

13.4 

14.0 

13.5 

11.7 

Equil. 3-pt. 

10.3 

11.0 

12.2 

11.0 

10.1 

Rhombic 4-pt . 

23.5 

25.2 

25.4 

25.2 

22.8 

Two-point; see Fig. E.5 

TQr . 

7.8 

18.5 

67.5 

3.5 

11.9 

TUr . 

67.7 

47.3 

88.2 

51.8 

13.8 

QrQr . 

46.2 

6.1 

33.3 

0.1 

3.6 

QrUr . 

92.4 

0.2 

94.5 

2.5 

0.1 

UrUr . 

75.8 

39.5 

58.0 

1.7 

3.8 


Pseudo-collapsed three-point; see Fig. E.5 


TTQr . 

69.3 

11.1 

33.4 

21.6 

11.1 

TTUr . 

5.3 

35.6 

70.1 

14.1 

13.7 

TQrQr . 

4.4 

11.0 

70.1 

24.5 

14.7 

TQrUr . 

39.6 

18.3 

98.5 

65.5 

58.0 

TUrUr . 

67.5 

2.6 

55.5 

1.2 

0.9 

QrQrQr . 

35.3 

5.3 

99.6 

44.0 

6.1 

QrQrUr . 

47.4 

0.2 

99.7 

41.7 

8.3 

QrUrUr . 

45.5 

62.4 

99.5 

14.6 

14.6 

UrUrUr . 

55.5 

60.5 

98.7 

75.7 

22.4 

Equilateral three-point; 

see Fig. E.5 





TTQr . 

1.6 

24.4 

50.0 

16.6 

18.7 

TTUr . 

67.4 

70.8 

81.5 

83.4 

68.3 

TQrQr . 

84.2 

14.5 

98.0 

56.3 

55.0 

TQrUr . 

20.3 

54.1 

95.6 

54.2 

30.0 

TUrUr . 

73.3 

2.2 

82.4 

22.2 

6.4 

QrQrQr . 

16.4 

53.1 

99.1 

6.5 

75.9 

QrQrUr . 

21.1 

80.5 

99.7 

93.4 

82.0 

QrUrUr . 

50.2 

46.1 

99.6 

59.7 

77.5 

UrUrUr . 

8.2 

34.0 

99.7 

88.1 

36.4 




Fig.E.6. Power spectra of the foreground cleaned CMB maps 
from FFP8 simulations. Top: TT power spectra evaluated us¬ 
ing the FFP8-UT74 mask. Bottom: EE power spectra evaluated 
using the FFP8-UP76 mask. Thick lines show the spectra of sig¬ 
nal plus noise estimated from the half-mission half-sum maps; 
thin lines show the noise levels from half-mission half-difference 
maps. The black line shows the input spectrum. 
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Fig.E.7. TT angular power spectra of residuals from the indi¬ 
cated FFP8 components in the Planck 2015 CMB maps, com¬ 
pared with the predicted signal from the best fit cosmology. 
“Other” is the sum of CO, free-free, thermal and kinetic SZ, 
spinning dust, and synchrotron emission. The horizontal axis is 
linear in £ . 


Fig. E.8. EE angular power spectra of residuals from the indi¬ 
cated FFP8 components in the Planck 2015 CMB maps, com¬ 
pared with the predicted signal from the best fit cosmology. 
“Other” is the sum of CO, free-free, thermal and kinetic SZ, 
spinning dust, far-infrared background, and radio and infrared 
unresolved sources. The horizontal axis is linear in . 
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